Looking into the hearts of Bok globules: 
MM and submm continuum images of isolated star-forming cores 
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ABSTRACT 

We present the results of a comprehensive infrared, submillimetre, and milUmetre continuum emission 
study of isolated low-mass star-forming cores in 32 Bok globules, with the aim to investigate the process 
of star formation in these regions. The submillimetre and millimetre dust continuum emission maps 
together with the spectral energy distributions are used to model and derive the physical properties of the 
star-forming cores, such as luminosities, sizes, masses, densities, etc. Comparisons with ground-based 
near-infrared and space-based mid and far-infrared images from Spitzer are used to reveal the stellar 
content of the Bok globules, association of embedded young stellar objects with the submm dust cores, 
and the evolutionary stages of the individual sources. Submm dust continuum emission was detected in 26 
out of the 32 globule cores observed. For 18 globules with detected (sub)mm cores we derive evolutionary 
stages and physical parameters of the embedded sources. We identify nine starless cores, most of which 
are presumably prestellar, nine Class protostars, and twelve Class I YSOs. Specific source properties like 
bolometric temperature, core size, and central densities are discussed as function of evolutionary stage. 
We find that at least two thirds (16 out of 24) of the star-forming globules studied here show evidence of 
forming multiple stars on scales between 1,000 and 50,000 AU. However, we also find that most of these 
small prototstar and star groups are comprised of sources with different evolutionary stages, suggesting a 
picture of slow and sequential star formation in isolated globules. 

Subject headings: dust, extinction — infrared: stars — ISM: clouds — stars: formation — submillimeter: ISM 



1. Introduction 

Different aspects of star formation can be stud- 
ied on different size scales and in different environ- 
ments. The large-scale distribution of star-forming re- 
gions and the relation between molecular cloud life 
cycles, galactic spiral density waves, and star forma- 
tion can be studied by observing nearby galaxies (e.g. 
Brouillet et aT][T998l [Scheepmaker et aLl[2009l l. The 
stellar initial mass function (IMF), which is needed 
to interpret these data, is usually derived from rich 



young stellar clusters in our own Galaxy (e.g. Kroupa 
|2002t [Chabrier|2003| l. Dense star-forming dark cloud 
complexes such as the Orion and Ophiuchus molecu- 
lar clouds are the places to study the relation between 
the molecular core mass spectrum (CMF) and the in- 



terstellar IMF (e.g..Motte et al.|1998| 


Nutter & Ward- 


Thompson|2007t|Simpson et al.|2008 


Goodwin et al. 
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Nearby and more isolated star-forming cores, such 
as Bok globules, are the best places to study in de- 
tail the initial properties of individual star-forming 
cores, their chemical evolution, kinematic structure, 
and the physics of their collapse and fragmentation 
(e.g., Clemens & Barvainis'1988';'Bourke et al. 1995a 
Launhai-dt & Henning 1997; Launhardt e t al.||1997' 
Henning & Launhardt 1998; Ward-Thompson et al. 
2007] [Stufz et al.||2008| l. Bok globules are small, 
simply-structured, relatively isolated, opaque molec- 
ular clouds that often contain only one or two star- 
forming core. They are often not completely iso- 
lated, but reside in the filamentary outskirts of larger 
dark cloud complexes ( Launhardt & Henning|1997 1, a 
fact that may tell something about their origin. With 
their size, mass, densities, etc., Bok globules resemble 
small clumps in larger molecular clouds (cf. Bergin 



& Tafalla"2007 ), only that they lack the surrounding 
cloud. Table [T] summarizes the average general prop- 
erties of typical Bok globules and their star-forming 
cores. 

Although they are the most simple star-forming 
molecular clouds, many globules deviate considerably 
from spherical geometry. They are often cometary or 
irregularly shaped. The dense star-forming cores are 
not always located at the center of the globule, but 
in cometary-shaped globules are often located closer 



'Dept. Experimental Physics, National University of Ireland 
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to the sharper rim at the "head" side (e.g., Launhardt 
1 1996 !). Similarly, pre-stellar cores, which are the ear- 



Uest stage of star formation (Ward-Thompson et al. 
1994 2QQ7\ often appear to be fragmentary and fil- 
amentary. However, the protostellar cores and en- 
velopes of the more evolved Class & I YSOs (| Lada| 



1987 Andre et al. 1993 1 are more spherically symmet- 
ric, which can be understood as a result of the gravi- 
tational collapse of the inner dense R ~ 5000 AU re- 
gion. Many of these isolated cores were found to be the 
sources of bipolar molecular outflows, indicating the 



presence of embedded protostars (e.g., |Yun & Clemens 
1994b| [Ward-Thompson eraL]|T995| [Bontemps et al. 
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In order to investigate the star-forming potential 
and evolutionary stages of Bok globules in the solar 
neighbourhood, we had surveyed a large number of 
globules for signs of star-forming cores, using as trac- 



ers, e.g., the mm dust continuum emission (Launhardt 



Tgge'; Launhardt & Hennin g] 1997[ [Henning & LauiT 
hardt 1998; Launhardt et al.|1998|l, NH3 (Bourke et al. 



1995b} or CS Une emission ( [Launhardt et al.|1998| l. 

In this paper we present a submillimetre and mil- 
limetre - hereafter (sub)mm - continuum study of 32 
Bok globules, which were identified from these previ- 
ous surveys as promising candidates for globules with 
currently ongoing star formation. The (sub)mm maps 
are complemented by deep near-infrared (NIR) images 
and NIR to mm spectral energy distributions (SEDs). 
In Sect.[2] we describe the observations and data re- 
duction. In Sect.[3] we present the NIR and (sub)mm 
images, SEDs, and results on multiplicity, physical pa- 
rameters, and evolutionary stages. In Sect.|4| we dis- 
cuss in particular the source properties as function of 
evolutionary stage and the results on multiplicity. In 
Sect.|5] we describe and discuss the individual glob- 
ules, and Sect.|6| summarizes the main results of this 
study. 

2. Observations and data reduction 

A total of 32 Bok globules were studied using a 
combination of NIR and (sub)mm continuum obser- 
vations. Coordinates, distances, and observation sum- 
mary of these sources are given in Table |2| Although 
many of the globules appear in different catalogues 
and thus often have several names assigned, we use 
throughout this paper the CB name and number in- 



troduced by [Clemens & Barvainis ( 1988 1 for northern 
sources, and the BHR name and number introduced by 
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|Bourke et al.| ( |1995al l for southern sources. 

2.1. (Sub)millimetre continuum observations 

Submillimetre and millimetre dust continuum ob- 
servations at 450 fj.m, 850 /im, and 1 .3 mm were per- 
formed during several observing runs between 1994 
and 2002 at the 15-m James Clerk Maxwell Telescop^ 
(JCMT) on Mauna Kea, Hawaii, at the 30-m Institut de 
Radio Astronomic Millimetrique (IRAM) telescope on 
Pico Veleta, Spain, and at the 15-m Swedish-ESO Sub- 
millimetre Telescope (SEST) on La Silla, Chile. 

Observations at 1.3 mm of the northern sources 
were carried out during four observing runs between 
April 1994 and March 1997 using the MPIfR 7- 
channel (1994-95) and 19-channel (1996-97) bolome- 
ter arrays ( [Kreysa et aLil998 ) at the IRAM 30-m tele- 
scope. The maps were obtained with the standard 



double-beam technique described by Emerson et al. 
1979| l. The telescope was scanned continuously at 



4"/sec in azimuth while the wobbling secondary, op- 
erating at 2 Hz with a chopper throw of 30" along the 
scan direction, generated the dual beams. Scan lengths 
and map sizes were adapted to source sizes to pro- 
vide sufficient basline for background determination. 
Individual scan maps of a source, which were later 
combined, were obtained at different hour angles in 
order to reduce scanning effects. 

Observations at 1.3 mm of southern sources were 
carried out in September and November 2002 at the 
Swedish/ESO Submillimetre Telescope (SEST) at La 
Silla in Chile, using the 37-channel bolometer array 
SIMBA (SEST Imaging Bolometer AiTay). Since 
SIMBA was used without a wobbling secondary mir- 
ror, several fast-scanning maps (at BO'Vsec) were ob- 
tained for each field at different hour angles to re- 
duce scanning effects. Source BHR71 was already 
observed in 1995 with the MPIFR single-channel 
bolometer at SEST ( |Bourke et aT][T997] l. Individual 
maps at different hour angles were obtained with the 
same double-beam mapping technique used at the 30- 
m telescope, but by employing the SEST focal plane 
chopper. 

The submillimetre observations were carried out 
during three observing runs between August 1997 and 
September 2001 using the Submillimetre Common 
User Bolometer Array (SCUBA) on the JCMT. The 
observations were conducted simultaneously at 850 



and 450//m using a 64-point jiggle pattern to create 
fully sampled maps of 2.3' diameter fields around the 
coordinate centers (for more details, see [Holland et al.| 
19991. Time-dependent variations in the sky emis- 



sion were removed by chopping the secondary mir- 
ror at 7.8 Hz with a chop throw of 120" in azimuth. 
During the 2000 and 2001 observations, we also ob- 
tained 850 jum polarization maps of six globule cores 
(see |Henning et al.||2001[ [Wolf et al.||2003 ). In addi- 
tion, the JCMT archive, operated by the Canadian As- 
tronomy Data Centre (CADC), was searched for addi- 
tional data for the sources observed at 1.3 mm. One 
of the sources (CB224) was observed using SCUBA 
in the scan-mapping mode. In this observing mode, 
the array is rastered across the sky to produce a rect- 
angular map several arcminutes in extent. Variations 
in sky emission are again removed by chopping the 
secondary mirror, though six chop configurations are 
used, chopping in 30, 44 and 68" in both RA and Dec 
Penness et al.|1998| ). 

During all observing runs, the atmospheric trans- 
mission was determined using the 'sky dip' method, 
which were performed every 2-3 hours. At the JCMT, 
we also used polynomial fits to the 225 GHz atmo- 
spheric opacity data which was measured by the Cal- 
tech Submillimeter Observatory (CSO) radiometer 
Average zenith optical depths for each run are listed in 
Table E 

Telescope pointing and focus were checked regu- 
larly by observing strong standard continuum point 
sources. Pointing was found to be repeatable within 
^3" at the IRAM 30-m telescope and the JCMT, and 
^5" at SEST. For flux calibration, maps of Uranus, 
Mars, Neptune, or different secondary calibrators 



(Sandell 1994) were obtained during each shift and 



at different elevations. Observing parameters for all 
runs are summarized in Table[3] 

2,2. Reduction and calibration of the (sub)millimetre 
data 

1.3 mm bolometer maps observed using the IRAM 
30-m and SEST telescopes were reduced and an- 
alyzed with the MOPSPl software package. The 
chopped dual-beam scan maps were flat-fielded, cor- 
rected for atmospheric extinction, de-spiked, base- 
line subtracted, and processed with an advanced sky 



'JCMT is operated by the JAC, Hawaii, on behalf of the UK STFC, 
The Netherlands OSR, and the Canadian NRC. 



^MOPSI (Mapping On-Off Pointing Skydip Infrared) is a software 
package for infrared, millimetre, and radio data reduction that has 
been developed and is updated by R. Zylka. 



3 



noise algorithm. The double-beam maps were then 
restored into single-beam maps using the Emerson- 
Klein-Haslam algorithm (Emerson et al. 1979| l. The 
individual maps were averaged with corresponding 
weights and converted into the equatorial coordinate 
system. The SIMBA fastscan maps were reduced ac- 
cording to the instructions in the SIMBA Observer's 



Handbook (2002; see also Chini et al. 2003 Appendix 
A). 

The SCUBA 450yum and 850//m jiggle maps were 
reduced with the SURF package (SCUBA User re- 
duction Facility, |Jenness & Lightfoot||1998| l. The in- 
dividual maps were flat-fielded, corrected for extinc- 
tion, processed with a sky-noise routine (Jenness et al. 



[1998 ), and combined. The SCUBA scan maps were 
processed in a similar manner, before removing the 
dual-beam function. Finally, the maps from the six dif- 
ferent chop configurations were combined in Fourier 
space, weighting each to minimize noise on spatial 



scales corresponding to the chop throw ( Jenness et al. 
[T998l|Hollandetal.|1999D . 

The data were calibrated using nightly observations 
of either primary planetary calibrators, or secondary 
calibrators when a planet was unavailable. Main beam 
flux calibration factors for all maps were computed 
by integrating over the main antenna beam in the 
planet maps, assuming planet brightness temperatures 
listed by [Orton et al.| ([T986| and |Griffin & Ortoiil 
( 1993)) and correcting for the error beams (see below). 
The flux densities of the secondary calibrators were 
taken from |Sandell| ( |1994| l. The statistical scatter of 
the 1.3 mm flux calibration factors derived from all 
IRAM 30-m and SEST data was found to be ^6%. 
Due to the uncertainty of the adopted planet temper- 
atures and uncertainties in the flux integration proce- 
dure, we estimate that the total calibration uncertainty 
is +15% at 1.3 mm, +20% at 850//m, and +30% at 
450/vm. Typical Icr noise levels (off-source) in the fi- 
nal maps are 5.5 + 1.3mJy/beam (IRAM 30-m) and 
20 ± 5 mJy/beam (SEST) at 1.3 mm, -20 mJy/beam at 
850 fim, and ~200 mJy/beam at 450 //m, respectively. 

For each instrument and observing run, the aver- 
age beam size full-width at half maximum (FWHM) 
was determined from the flux calibration maps, mostly 
on Uranus, taking into account the intrinsic size of the 
planet disk (see Table|3]l. The secondary calibrators are 
less well characterized, and so on nights when plane- 
tary beam-maps were unavailable, the average beam 
FWHM for the entire data-set is quoted in Table |3] 
Error beams at 1.3 mm and 850 jum were found to be 



negligible, but at 450 yum a significant error beam was 
found to contribute up to ~50% of the flux density. 
Note that the main beam flux calibration accounts for 
this (to first order). Flux density measurements in the 
final maps were done using the GILDASj^software and 
are are described in Sect.l3.1l 



For weak and extended sources, the integrated flux 
density and the derived source size also depend on 
the map size, the curvature of the extended flux den- 
sity distribution, and the method of baseline subtrac- 
tion, since the sources are always observed against a 
more extended emission background originating from 
the parental molecular cloud. 

2.3. Near-infrared observations 

We use the NIR images only for overlays with the 
(sub)mm maps, to identify which core is associated 
with a star or NIR nebula, and for discussion of evo- 
lutionary stages. NIR photometry is not be described 
in this paper, but we use the fluxes for the SEDs where 
relevant. To this end, NIR observations of the globules 
listed in Table [2] were performed during five observ- 
ing runs between December 1993 and October 2003 
at the 3.5m telescope on Calar Alto, Spain and at the 
MPG/ESO 2.2-m telescope on La SiUa, Chile. Observ- 
ing parameters for all runs are compiled in Table |4] 
Images were obtained in standard dithering mode with 
4 to 5 positions per field. Flat field calibrations were 
obtained each night via dome flats using difference 
frames of integrations with a halogen lamp on and off. 



Sets of photometric standard stars from lists by EUas 



let al.| ( |1982| i and |Hawarden et al.| ( |200T] i were observed 
during each night, but the calibration for most sources 
was done by using 2MASS catalogue fluxes for stars 
in the fields. 

The data were reduced using a combination of the 
IRAlj^ "Experimental Deep Infrared Mosaicking Soft- 
ware" (XDIMSUnQ and Starhnk Project packages. 
For photometric calibration, zero-points were derived 
from observed standard stars of each night. Then, 
magnitudes for individual objects were derived using 



^http://www.iram.fr/IRAMFR/GILDAS 

''iRAF is distributed by the NOAO, which are operated by the Asso- 
ciation of Universities for Research in Astronomy, Inc., under coop- 
erative agreement with the National Science Foundation. 

^XDIMSUM is a variant of the DIMSUM package (V3: August 19, 
2002) developed firstly by Eisenhardt et aT7|p996| 
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GAIA0with the embedded PHOTOlvQ package. The 
average 3 cr point source detection Hmit in all three 
NIR bands is ~ 19 mag. Only sources that are associ- 
ated with a dense globule core (i.e. a (sub)mm source) 
are considered in this paper. 

For astrometric calibration of mosaicked frames we 
used the Guide Star Catalogue i:0and the SIMBAE^ 
The resulting typical positional accuracy is better (X'S 
(1 cr rms). 

2.4. Spitzer data 

All sources listed in Table|2] were observed by 
Spitzer within different programs, most with IRAC 
and MIPS, except CB 17 SW, BHR58, BHR 137, and 
CB 224. The relevant data were obtained from the 
Spitzer Science Center (SCC) Archive. In particular, 
we retrieved, where available, IRAC maps at 3.6, 4.5, 
5.8, and 8;um and MIPS scan maps at 24, 70, and 
160 jum. The data were processed by the SCC using 
their standard pipeline to produce Post-Basic Cali- 
bration Data (P-BCD) images. Photometry was done 
in IDL by applying aperture photometry, where the 
aperture corrections available at the instruments web- 
page were used. The conversion from P-BCD units 
(MJysr"') to flux density was achieved by applying 
Quick Point Source Photometry as described online. 
The overall flux calibration and flux measurement ac- 
curacy is estimated by the SSC to be accurate to within 
20%. 

3. Results 

3.1. (Sub)millimetre continuum maps 

In Figs. [T| through 24 we present the (sub)mm 
maps, together with NIR K-band images for each of 
the sources listed in Table |2] (sources with only upper 
hmits ai-e not shown). Digital Sky-Survey (DSS2|^ 



*GAIA is provided by Starlink, and is a derivative of tlie Skycat cata- 
logue and image display tool, developed as part of tlie VLT project at 
ESO. Skycat is free software under the terms of tlie GNU copyright. 
^The PHOTOM photometry package is provided by Starlink 
^The Guide Star Catalogue II is a joint project of the STSI and the 
Osservatorio Astronomico di Torino. STSI is operated by the As- 
sociation of Universities for Reseai'ch in Astronomy, for the NASA 
under contract NAS5-26555. 
'SIMBAD is an online database, operated at CDS, Strasbourg, France 
'"The Digitized Sky Survey was produced at the Space Telescope Sci- 
ence Institute under US Government grant NAG W-2166. The im- 
ages of these surveys are based on the photographic data obtained 
using the Oschin Schmidt Telescope on Palomar Mountain and the 



optical images of the globules are also shown, together 
with the submm contours (850yL(m or 1.3 mm) and the 
boundaries of the submm maps, to illustrate their over- 
all morphology and the relative location of the dense 
cores within the globules. Note that in the case of 
SCUBA jiggle maps the entire mapped region is indi- 
cated, while in the case of IRAM (dual-beam) or SEST 
scan maps only the area from which the mm emission 
was restored is marked. In addition, SEDs, includ- 
ing fluxes from optical wavelengths up to 3 mm, are 
shown for those sources were we could compile suf- 
ficient data and associate them unambiguously to one 
source. 

(Sub)mm flux densities are summarized in Table [5] 
The center coordinates for each dust core were derived 
from gaussian fits to the 1.3 mm intensity maps af- 
ter decomposition into extended envelope and individ- 
ual compact components. For some sources (marked 
in Table |5]l, the peak position was derived from the 
850//m maps. The peak intensity ll'^'^^ of each com- 
pact component is given in Columns 3, 5 and 7 of Ta- 
ble[5]in units of Jy/beam. The beam size is also quoted 
here in parentheses. Note that the maps of some ex- 
tended sources (e.g., CB 246) were convolved with a 
larger gaussian beam to increase the SNR. We also 
quote the total integrated flux density, 5™, for each 
source. This was derived by integrating within a poly- 
gon enclosing the (closed) 2 cr contour. If there is one 
or more compact sources within a larger envelope, we 
also give the integrated flux densities of the compact 
sources, derived after fitting and subtracting the ex- 
tended envelope from the maps. We do not list indi- 
vidual error bars for the flux densities, because uncer- 
tainties are mostly dominated by the calibration uncer- 
tainties (see Sect. |2.2| i. However, for extended sources, 
the uncertainty is dominated by the size of the area 
from which emission could be recovered above the 2- 
cr-level, which in turn depends on the curvature of the 
emission background, the ratios of the angular sizes of 
the beam, the source, and the chop throw, as well as on 
the SNR in the map. In general, the relative fraction of 
recovered to total actual flux increases from faint ex- 
tended to strong compact, and from nearby to distant 
sources. 

3.2. Association of sub(mm) and infrared sources 

When compiling SEDs and attempting to model 
and interpret the observational results, it is important 

UK Schmidt Telescope. 
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to reveal which of the sources and flux densities de- 
tected at different wavelengths, with different angular 
resolution, apertures, and astrometric precision, are ac- 
tually physically associated and which ones have dif- 
ferent origins. For this reason, we overlay the (sub)mm 
maps on the optical (DSS) and NIR images, show the 
IRAS point source position error ellipses where appli- 
cable, and also inspect the Spitzer maps where avail- 
able (see \2A\ . 

We found the following associations between 
NIR/MIR sources and mm cores: No NIR sources 
brighter than 0.02 mJy (s; 19 mag) at 2.2 fim or ^ 1 mJy 
at 4.5 /im are detected in CB 17-SMMl and 2, CB 130- 
SMM2, CB232-SMM, CB243-SMM1 and 2, and 
CB 246. However, since local dust extinction, source 
morphology, and projection effects (e.g., the ahgnment 
between the line of sight and outflow cavities) play a 
major role at these early evolutionary stages, this de- 
tection Umit cannot be directly converted into an upper 
Umit for the mass or luminosity of possible embedded 
sources. No K-band emission, but faint NIR sources 
with steeply rising NIR SEDs peaking at a; 4 ... 6 /im 
have been found in BHR55, BHR86, CB68, and 
CB 199. Extended NIR nebulosities, in several cases 
NIR jets, but no star-like objects were detected in 
CB 26, BHR 1 2-SMM 1 and SMM2, BHR 36, BHR 7 1 , 
CB 130-SMMl, CB 224-SMMl, CB 230, and CB 244- 
SMMl. A single red star or compact star-like NIR 
source, in some cases with extended nebulosity, was 
detected at the mm peak position in CB 6, CB 188, and 
CB 224-SMM2. Several NIR sources or clusters, re- 
spectively, are associated with the (sub)mm cores in 
CB 34, CB 54, CB 205, and CB 240. The situation is 
somewhat inconclusive in CB 58 and BHR 137, where 
it is not clear if the NIR sources are actually associated 
with the mm core or not. 

The following (sub)mm cores have very nearby, but 
most likely not directly associated NIR/MIR sources 
and the flux measurements could not always be dis- 
entangled: CB 17 has an IRAS point source and a 
5p/fzer-detected low-luminosity IR point source that 
is clearly not associated with the mm core (Fig.|2|i. 
CB 130-SMMl is a strong and slightly extended 
submm source with a small NIR nebula at the mm peak 
position. A very nearby (14") strong and pointlike 
IR source (IRS) could not be completely disentagled 
from SMMl, but does not seem to exhibit detectable 
mm emission on its own. CB232 and CB243 both 
have bright reddened stars within the area of the mm 
emission, but clearly offset and not associated with the 



main mm peak. 

3.3. Multiplicity 

Although the (sub)mm maps are not sensitive to 
spatial scales that cover the peak of the binary star 
separation distribution (20-lOOAU; Duquennoy & 
|Mayor|199"T] l, we find that at least 16 out of 24 glob- 
ules with detected (sub)mm emission show signatures 
of multiplicity at scales larger than a!l,OOOAU, out 
to 80,000 AU (0.4 pc; see Table|6]l. This includes 
prestellar and protostellar double cores resolved in the 
mm dust continuum (e.g., CB 17, CB 130, CB246, 
BHR 12), unresolved protostellar cores where NIR or 
MIR observations have revealed a secondary embed- 
ded source (e.g., CB230, BHR71), prestellar or pro- 
tostellar cores where NIR/MIR observations have re- 
vealed sources that are significantly offset from the 
(sub)mm peak and that are most likley not embedded 
in these cores (e.g., CB 17, CB 130, CB 232, CB 243; 



see Fig. 28 1, as well as large dust cores that are associ- 
ated with clusters of NIR sources, presumably young 
stars that formed from the same globule (e.g., CB34, 
CB 54, CB 205). Nine of these 24 globules were re- 
solved to have two or more cores in the (sub)mm 
dust emission maps (see Table[6]l. However, only 
three of these double cores have projected separations 
<10,000 AU (CB 17, BHR 12, and CB 130). Separa- 
tions > 10,000 AU are probably not relevant for the for- 
mation of bound stellar systems. 

Altogether, we find that at least two thirds (16 out 
of 24) of the star-forming globules studied here show 
evidence of forming multiple stars, either in multiple 
star-forming cores, wide embedded binaries, or small 
star clusters. The fraction of closer binaries formed 
from unresolved (sub)mm cores might be higher, but 



remains unknown from this study (cf. Maury et al. 
|20T0| . 

3.4. Spectral Energy Distributions 

Spectral Energy Distributions (SEDs) were com- 
piled from aperture photometry in the (sub)mm maps 
(Table|5]l and NIR images (Figures [T 24 1 presented 
in this work (Sects. |TT] and |12|, the 2MASS cat- 
alogue, the IRAS point source catalogue, pipeline- 
calibrated Spitzer IRAC and MIPS maps, and other 
publicly available databases and literature papers. We 
do not explicitly list all flux values which we have 
compiled, but show in Figs.[T]through 23 the SEDs for 
those sources where we are confident that the measure- 
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merits represent one source or, in some cases, physi- 
cally related, mostly unresolved double sources (e.g., 
CB 230). Note that we actually show the energy den- 
sity (vSv) and not the flux density (Sy)- We do not 
show SEDs for sources where we either do not have 
enough data points (e.g., BHR 137, CB 246) or where 
the (sub)mm emission is associated with multiple NIR 
sources (e.g., CB 34, CB 54, CB 58, CB 205, CB 240). 
For the other sources, we indicate in the figure captions 
and in TablejTjwhether we use the total (sub)mm fluxes 
or only one compact component. Non-detections are 
indicated by down-arrows as 3 cr upper limits. Flux 
values that derive from apertures smaller than the ex- 
pected emission region (e.g., single-point measure- 
ments with the JCMT-UKT14 bolometer) or interfer- 
ometric measurements, which are not sensitive for ex- 
tended emission, are indicated by up-aiTows as lower 
limits. A comprehensive overview of all SEDs, or- 



dered by evolutionary stage (see Sect. 3.5 1 is shown in 



Fig. 25 



In order to derive luminosities and bolometric tem- 
peratures, we first interpolated and then integrated the 
SEDs, always assuming spherical symmetry. In the 
case of, e.g., bipolar NIR reflection nebulae, this may 
lead to either an under-estimate of the NIR luminos- 
ity (e.g., CB 26, Fig. |3]l or an overestimate (e.g., 
CB230, Fig. 19 1, depending on the inclination with 



respect to the line-of-sight. However, all sources stud- 
ied in this paper emit the majority of their total en- 
ergy at much longer wavelengths, where the radiation 
is reprocessed by dust, so that the general results de- 
rived here are not significantly affected by the sim- 
plifying assumption of spherical symmetry. Interpo- 
lation between the measured flux densities was done 
by ax^ grey-body fit to all points at A > 100//m and 
by simple logarithmic interpolation between all points 
at A < lOOyum. Free parameters in the grey-body 
fits were the dust temperature, source size, and op- 
tical depth. The (sub)mm spectral index of the dust 
emissivity was fixed to p 



3.61. Fit- 



1.8 (see Sect. 

ted dust temperatures were in the range 15 K (CB 17- 
SMM) to 53 K (CB 224-SMM2), with mean values of 
33 K for Class sources and 42 K for Class I sources 
(see Sect.[33]l. 

Note, however, that we use these fits only to de- 
rive an approximate interpolation for integrating the 
SEDs. In particular, we do not use the dust temper- 
atures or source size from these grey-body fits because 
the derived effective temperatures are usually dom- 
inated by small contributions from either embedded 



heating sources or externally heated very small grains 
in the envelopes. Thus, they usually do not repre- 
sent the mass-averaged dust temperature. Instead, we 
use fixed mass-averaged dust temperatures as function 



of evolutionary stage (see Sect. 3.6 1 and derive source 



sizes directly from the (sub)mm maps. 

3.5. Evolutionary stages 

We use a combination of different diagnostics to 
derive the evolutionary stages of the globule cores. 
We do not attempt to classify those cores of dis- 
tant globules that are associated with clusters of NIR 
sources and for which the NIR, FIR (IRAS), and 
(sub)mm emission could not be uniquely associated 
(e.g., CB 34, CB 54, CB 58, CB 205, CB 240). Sources 
for which we have barely more information than a mm 
flux (e.g., BHR 58, BHR 137) ai-e also not considered 
here. For the other sources, we analyze the following 
tracers of evolutionary stage: 

1. The ratio of submm to bolometric luminos- 
ity, Lsmm / Lboh where Lsmm is measured long- 
ward of 350 fim, which roughly reflects the ratio 
of envelope to stellar mass, Menv / Mt ( [Andre 



et al.|1993| l. [Andre et al.] ( |2000| l propose to use 
^smm / L\,o\ > 0.5% as a threshold for Class 
protostars, reflecting Menv / M» > 1 . 



2. 



The bolometric temperature, Tboi (Myers & 
Ladd 1993[ ), which is defined by the tempera- 



ture of a blackbody having the same mean fre- 
quency (ratio of the first and the zeroth moment 
of the SED) as the observed continuum spec- 
Chen et al.[ ( |T995] l and [Myers et"aL| ( [TW8| ) 



trum. 



derive the boundary between Class and Class I 
YSOs to be at Tboi ~ 70 K. Note that we use the 



interpolated SEDs (see Sect. 3.4 1 with flux val- 
ues at all wavelengths to compute Tboi and lumi- 
nosities, rather than just the actual data points. 
The latter approach could lead to significant bi- 



ases if the SED is not well-sampled (cf. Enoch 
[et al.|2009l l. 

3. The IR spectral index, cir = -d log(v Fy) / c/ log v 
( [Lada| 1987) 1, which characterizes the IR excess 
with respect to stellar blackbodies. Where suf- 
ficient data are available, we calculate ajR be- 
tween A - 2.2 and 24 yum, thus avoiding to car- 
acterize only the NIR/MIR nebulosity "bump" 
in the Class SEDs (see below), the possible 
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influence of the silicon absorption, and the low- 
resolution IRAS measurements at 12 and 25 fim, 
which are prone to confusion problems. Quanti- 



tatively, we follow Andre & Montmerle ( 19941 
and use aiR = as boundary between Class! 
and II. No well-defined boundary exists for 
Class sources since only Spitzer lead to the 
detection of the NIR/MIR emission in many of 
them. 

Along with these three quantitative classification 
criteria, we evaluate the morphology of the NIR emis- 
sion and the presence of molecular outflows. The first 
one can be either direct star light, indicating that the 
circumstellar emission originates either from an op- 
tically thin envelope or a phase-on disk, or diffuse 
nebulosity, indicating, e.g., scattered light from out- 
flow cavities in an otherwise optically thick envelope. 
The presence of a molecular outflow is an indicator 
of an embedded source that distinguishes protostars 
from prestellar cores. Two notes of caution are: (i) the 
non-detection of an outflow does not prove the absence 
of an embedded protostar, in particular as we have no 
homogeneous high-sensitivity, high-angular resolution 
outflow survey at hand, and (ii) the driving source of an 
outflow is not necessarily identical with the (sub)mm 
source - this must be verified in each individual case. 

The applicabilty and drawbacks of the classification 
criteria mentioned above has been discussed in detail 
by [Enoch et al.|p009| l. An analysis of our data con- 
firms the finding of these authors in that does not 
correlate well with Tboi and embeddedness for Class 
and embedded Class I sources. The deep NIR images, 
in particular those from Spitzer, show that most of even 
the youngest Class sources already have a small NIR 
nebulosity associated. This can also be clearly seen 



in the compilation of SEDs shown in Fig. 25 We at- 
tribute this to the fact that the amount of NIR/MIR flux 
seen from such an embedded object largely depends on 
the alignment between the line of sight and the outflow 
cavity, allowing or not a direct view onto the hot dust 
in the inner protostellar core and/or forward-scattered 
light in the outflow cavity. Therefore, we only list the 
aiR values in Table]?] but do not use them for the clas- 
sification of evolutionary stages. 

We find that all sources, for which we could com- 
pile meaningful SEDs, have L^rnm I L\kA ^ 0.8% (Ta- 
ble[7)l, i.e., according to the criterion proposed by 



ilar properties of the SEDs, appears to be a more 
reliable indicator of the evolutionary stage than the 
L^mm I ^'boi ratio, as it is less sensitive to the amount of 
extended envelope emission recovered in the submm 
maps (see Sect. 4.1 1. We therefore use the bolometric 



temperature as the main discriminator between Class 
protostars and Class I YSOs, along with the other sup- 



porting indicators mentioned above (cf. jEnoch et al. 
[20091 ). 

The derived evolutionary stages, along with other 
source properties, for 18 globules are summarized in 
Table]?] Figure]26] shows the Lsmm / ^boi ratio vs. Tboi 
for these sources. In total, we identified nine star- 
less cores, presumably prestellar, in six globules, nine 
Class protostars in nine globules, and twelve Class I 
YSOs in eleven globules. Interestingly, seven glob- 
ules contain co-existing sources at different evolution- 
ary stages (see discussion in Sect. [4.2| i. Six globules 
were not detected at (sub)mm wavelengths, and for the 
submm cores in eight globules we could not or did not 
attempt to derive an evolutionary stage because we ei- 
ther had insufficient data, or the available data indi- 
cated that several sources with possibly different evo- 
lutionary stages contribute to the partially unresolved 
flux measurements. 

3.6. Physical parameters of the star-forming cores 

The following physical source quantities are de- 
rived from the (sub)mm continuum maps and the 
SEDs: bolometric and submm luminosities (Sect. [3.5l l, 
source sizes, masses, beam-averaged peak column 
densities, and source-averaged volume densities. 

FWHM mean source sizes ( Vfll?) are simply de- 
rived by measuring the angular diameters of the 50% 
contours in the 850 jum or 1.3 mm maps (after decom- 
position into envelope and compact core, where appli- 
cable), deconvolving them with beam size, and apply- 
ing the distances listed in Table]2] The results are listed 
in Table]?] 

To calculate masses, densities, and column densi- 
ties, we adopt mass-averaged mean dust temperatures 
that depend on the evolutionary stage of the source 



(Sec t.|3.5|l . Quantitatively, we follow [Motte & Andr6 



Andre et al. ( 2000| they would all be Class sources. 
The bolometric temperature, although reflecting sim- 



( |200l| , lShirley et al.[ ( [2052l i, and [Young et al.[ ( |2003| l 

who have derived mean mass-averaged isothermal dust 
temperatures from radiative transfer models for sam- 
ples of starless cores. Class protostars, and envelopes 
of Class I YSOs. We use here {T^) = lOK for star- 
less (prestellar) cores, {Ti) = 15K for Class pro- 
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tostars, and (To) - 20K for the envelopes of Class I 
YSOs. The first dust temperature map directly derived 
from new Herschel data of the Bok globule CB244 
with its two embedded cores confirms the validity of 
these temperature assumptions fS tutz et al.|2010) l. The 
adopted temperatures and resulting masses for each 
source are listed in Table|7] 

Assuming optically thin isothermal dust emission, 
the hydrogen gas mass Mh = M(H) + 2 M(H2) is cal- 
culated from the (sub)mm flux density by 



(1) 



Here, v is the observing frequency, Sv the observed 
flux density from the region of interest, D is the 
distance to the source, K^iv) is the dust opacity, 
Byiv, Ti) the Planck function, the dust tempera- 
ture, and Mh/Mj the hydrogen-to-dust mass ratio. 
Throughout this paper, we use a hydrogen-to-dust 
mass ratio of Mh/Mj - 110 and a dust opacity of 
ka{A) - /Cd(1. 3mm)x (/l/1.3mm)~''^ w ith 1.3mm) = 
O.Scm^g ' (Ossenkopf & Henning 1994. uncoagu- 



lated MRN dust grains with thin ice mantles), a that 
is value typical for dense and cold molecular cloud 
cores and intermediate between the opacities usually 
adopted for unprocessed interstellar grains in the dif- 
fuse gas (e.g., |Draine & Lee| 1984 1) and for processed 
and coagulated grains in protoplanetary disks (e.g., 
IBeckwith et al.|2000| i. 

We derive total Hydrogen masses of 1-20 M© for 
prestellar cores, 0.5 - lOM© for Class protostars, and 
0.11 -IMq for Class! YSOs (not counting the ex- 
tended envelopes; see Table[7]i. The total gas mass, 
accounting for helium and heavy elements, amounts 
then to Mg ^ 1.36 Mh. 

The beam-averaged peak Hydrogen column density 
Nyi - N{\i) + INiHi) was derived from the observed 
peak flux densities by 



A^H = 



1 /Mh\ 

— h# ' (2) 



Ki{v)^lsBy{v,Ti) mn \Md 



where ly is the peak flux density (Table|5]l, Q.b = 
1.133 0^ the beam soUd angle, flg the beam FWHM 
size, and is the the proton mass. The resulting 
beam-averaged peak column densities are also listed 
in TableEl 

Assuming that the sources have the same extent 
along the line of sight as the mean size perpendicular 
to the line of sight, the source-averaged volume density 



of hydrogen atoms, hh = "(H) + 2 n(H2), was calcu- 
lated by 

"H = 77 , (3) 

mn V 

with V ~ 7r/6t/j being the volume of the source, and 
ds the FWHM size listed in Table |7] The resulting 
densities are listed in Table|7] Note that these are 
source-averaged densities and the local peak densities 
in sources with density gradients can be orders of mag- 
nitude higher. 

In an earlier paper, we had also derived radial den- 
sity profiles for most of the globule cores by fitting ray- 
traced and beam-convolved model cores to the radially 
averaged (sub)mm data at all three wavelegths ( jLaun-] 
|hardt|2005) . The extended dataset presented in this pa- 
per would certainly permit to derive better-constrained 
density profiles for many globules cores. However, as 
in this previous paper, such an analysis would be ham- 
pered by the fact that we have to assume temperature 
profiles, since the (sub)mm data provide only weak 
constraints. In particular for the prestellar cores, which 
do not have a central heating source and may have a 
positive radial temperature gradient, this can introduce 
significant biases and may lead to an underestimation 
of the density gradient. Since we expect to obtain soon 
Herschel data for many of the sources presented here, 
which will provide direct constraints on the tempera- 
ture profiles, we postpone the derivation and analysis 
of the radial density profiles to a succeeding paper. 

4. Discussion 

4.1. Source properties as function of evolutionary 
stage 

Being relatively nearby and isolated, the cores of 
Bok globules are the best-suited places to observe the 
properties and evolution of individual low-mass star- 
forming cores. However, this isolation, which makes 
them such good targets, may also alter the core proper- 
ties and evolutionary process as compared to the more 
representative (for general star formation), but harder 
to observe, star-forming cores in clustered regions. 
Some of the most obvious differences in environmen- 
tal conditions are the interstellar radiation field, which 
is not shielded by the surrounding cloud, the miss- 
ing external pressure due to the lack of inter-core gas, 
and the lack of turbulence input from the surrounding 
cloud. 

The most obvious changes in core properties during 
the evolution from prestellar cores, through collapsing 
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protostars, towards Class I YSOs are expected to be 
seen in the core size, central core density, ratio of core 
mass to envelope mass, and contribution of the central 
(proto)star to the source luminosity at shorter wave- 
lengths. The latter two effects have been used in the 
form of parameters Lsmm / ^boi and T\,oi as main cri- 
teria to identify the evolutionary stage of the sources 
(Sect.|33]l. 

In Figure 26 we show a plot of Lsmm/iboi vs. Tboi 



sources as "Class I" in Table]?] despite the fact that 
they have aiR < 0. 



for the sources compiled in Table]?] We find that all 
sources with Tboi < ?0 K, the empirical boundary be- 
tween Class protostars and Class 1 YSOs proposed 
by [Myers et"aL](]T998]l, have L^mm I i-boi > 2%. Within 



the Class phase, there seems to be no significant cor- 
relation between the two parameters and no obvious 
evolution from lower to higher values of Tboi- This 
may reflect the fact that during this early phase, the 
optically thick envelope dominates over the still very 
low-mass and low-luminosity central protostar. How- 
ever, there is a relatively large spread in Lsmm / Lboi ra- 
tios, probably reflecting the fact that during this phase 
the mass of the infalling envelope is relatively quickly 
reduced. 

The more evolved (Class 1) YSOs in our sam- 
ple have bolometric temperatures between lOOK and 
350 K and Lsmm/Lboi ~ 0.8 . . . 3.5% (see Fig.]26]). All 



sources with NIR nebulosities but no visible stars 
have Lsmm/Lboi ^ 1.6% and extended (resolved) 
(sub)mm emission. These sources are probably less- 
evolved, envelope-dominated Class 1 YSOs. The three 
sources with U^^/Uoi < 1.3% (CB 188, CB224- 
SMM2, CB232-1RS1) all have visible stars and com- 
pact (sub)mm emission arising presumably from cir- 
cumstellar disks. These sources probably represent 
the more evolved, disk-dominated Class 1 YSOs. We 
take these findings as tentative indication that the 
Lsmm / Lboi ratio, as indicator of envelope dispersal, 
may better trace the evolution within the Class and 
Class 1 phases than Tboi. 

Two of the Class 1 sources in our sample have 
double-peaked SEDs with NIR nebulosities and aiR < 
(CB26, CB230), which would indicate Class 11. 
However, the IR emission seems completely originat- 
ing from scattered light and the IR fluxes and morphol- 
ogy seem dictated by projection effects. For CB 26, 
a massive edge-on protoplanetary disk with a rem- 



nant envelope, Launhardt & Sargent (2001 1 derive an 
age of 1 Myr. CB 230, with a more massive enve- 
lope and less-developed disk(s) ( |Launhardt]|2001) l, is 
most likely younger. We therefore classify these two 



To verify these relations between Lboi, Lsmm/^^i 



boh 



and evolutionary stages, we have compared our sam- 
ple of globule cores to the sample of Class proto- 
stars from Andre et al. {ZQQQ) . Plotting these sources 



into the same diagramme (not shown here) shows a 
different picture. While none of the globules sources 
have both Tboi < ?0K and Lsmm/Lboi < 2%, numer- 
ous protostars from the Andre sample do populate this 
region. The two different source samples have similar 
distance distributions, and in each case, both nearby 
and more distant sources are distributed evenly around 
the Lsmm/i'boi VS. Tboi plot. Therefore there is no ev- 
idence to suggest that limited resolution or sensitivity 
for either of the samples is causing the discrepancy. 
Similarly, there is no evidence that any particular as- 
pect of the morphology of a Bok globule is causing 
this discrepancy between the two samples. As none 
of the above characteristics can possibly cause a sys- 
tematic shift of the Bok globules in the Lsmm/Lboi vs. 
Lboi diagramme, we conclude that the absence of glob- 
ule cores with both low Tboi and low Lsmm /Lboi -ratios 
is a real effect, caused by the systematic environmental 
differences between the cores of Bok globules and pro- 
tostars located in larger molecular clouds mentioned at 
the beginning of this section. 

The outer boundary of a protostellar core that is em- 
bedded in a larger molecular cloud is generally judged 
to be where it's envelope merges with the surrounding 
cloud material, or the envelope of an adjacent object. 
In contrast, the column density profile of a globule core 
continuously decreases towards the "edge" of the small 
cloud, making it more difficult to distinguish between 
the core and the surrounding cloud, both physically 
and observationally, as the envelope (cloud) emission 
is usually not chopped away in the small globules. This 
has the effect of making the detected envelopes larger 
and elevating the submm luminosity as compared to 
dense cores embedded in larger clouds. It is therefore 
possible that the cutoff in Lsmm/Lboh used to differ- 
entiate between Class protostars and Class 1 YSOs 
should be raised from 0.5% in the case of protostars 
in larger molecular clouds, to » 2% in the case of iso- 
lated Bok globules, to account for their larger enve- 
lope of cold dust. The value of 2% is chosen as this 
is the lowest value of Lsmm/Lboi in the observed sam- 
ple of globule cores with Lboi < ?0K. Although we 
cannot present a modified quantitative physical model 
(cf. |Andre et al.|1993| l, we fliink that this different cut- 
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off value for Lsmm/Lhoi still represents the point where 
the mass of the infalling protostellar envelope is equal 
to that of the accreting protostar, but it accounts for 
the contribution from the extended envelope which can 
no longer be clearly distinguished from the protostellar 
core. 

We conclude that Tboi appears to be the most re- 
liable parameter to discriminate between Class pro- 
tostars and Class 1 YSOs, with no hint for a system- 
atic difference between sources in isolated globules 
and those embedded in larger molecular clouds. On 
the other hand, the evolution during the Class phase 
and during the Class I phase may probably be better 
traced by the Lsmm / i'boi ratio. Despite the fact that the 
two parameters are not completely independent, this 
can be understood as a consequence of Tboi being more 
sensitive to a flux contribution at shorter (NIR) wave- 
length from a central heating source, while Lsmm / i^boi 
is more sensitive to the size and mass of the cold enve- 
lope. 

We have also determined the mean FWHM core 
sizes and the mean source-averaged volume densities 
(see Table|7]i of the three evolutionary subgroups. We 
find that the mean core sizes decrease from a;8,000 AU 
for prestellar cores, to =s3,000AU for Class proto- 
stars, to <2,000AU for Class 1 YSOs (Fig.W}. We 



also find that the Class 1 YSOs usually have extended 
envelopes with a mean diameter of 20,000 AU. The 
source sizes of the Class protostars might be some- 
what underestimated since we did not decompose the 
sources into central point source (the warm protostar) 
and envelope before measuring the FWHM size. The 
source-averaged volume densities, rtu, increase from 
1x10^ cm"-' for prestellar cores, to 7 x 10^ cm"^ for 
Class protostars. At these high densities, one can 
in general assume that gas and dust temperatures are 
coupled. For the cores of the Class 1 YSOs, which are 
mostly unresolved and presumably contain significant 
flux contrinutions from circumstellar disks, we derive 
a lower limit to the mean density of > 8 x 10^ cm"^. 
The extended remnant envelopes of the Class I YSOs 
have a mean density of 7 x lO'* cm"-'. 

4.2. Multiple and sequential star formation in 
globules 

One interesting aspect of the combined results for 
multiplicity and evolutionary stages is that the large 
majority of the globules with signs of multiple star 
formation shows evidence of non-coeval, possibly 



sequential star formation. This includes neighbour- 
ing mm sources with obviously different evolutionary 
stages (e.g., BHR12, CB 130, CB224), prestellar or 
protostellar cores with nearby IR sources, presumably 
more evolved protostars or Class 1 YSOs (e.g., CB 17, 
CB 130, CB 232, CB 243), as well as NIR star clusters 
next to large (sub)mm cores with the potential to form 
more stars (e.g., CB34, CB54, BHR137) (see Ta- 
ble|7]i. In only three globules we find (approximately) 
coeval pairs, ranging from multiple pre-stellar cores in 
CB 246, embedded Class protostars in BHR71, to a 
pair of embedded Class 1 YSOs in CB 230. 

The most interesting sources in this respect are 



shown in Fig. 28 CB 17 contains two prestellar cores 
and one low-luminosity Class 1 YSO within 5 000 AU 
(projected separation). CB 130 contains three sources 
within 9 000 AU, all with different evolutionary stages. 
A Class protostellar core is flanked by a prestellar 
core at 6,000 AU to the west, and by a Class 1 or II 
YSO at 3,000 AU to the east. CB 232 and CB 243 both 
contain a prestellar core with a Class 1 YSO within 
7000AU. 

This widespread non-coevality suggests that the 
multiple sources observed here have in general not 
formed by gravitational fragmentation of collapsing 
IBonneU & Batel '2006), 



protostellar cores (e.g.. 



but 



rather by initial turbulent fragmentation prior to the 
collapse phase (e.g., |Padoan & Nordluiid||2002| l, fol- 
lowed by a relatively independent evolution of the in- 
dividual fragments. Whether or not the first forming 
star can trigger collapse of the other fragment(s) and 
thus lead to small-scale sequential star formation can- 
not be directly answered by this study. The lack of 
close coeval pairs can be understood as an observa- 
tional bias as our study is in general not sensitve to 
size scales <1000 AU, which would be more relevant 
to the first mechanism. Consequently, we cannot draw 
direct conclusions on the frequency of bound binary 
stars formed in isolated Bok globues (cf. |Maury et al.| 
|20T0| . 

These findings also call for special attention when 
compiling SEDs and attempting to derive source prop- 
erties from flux measurements with insufficient angu- 
lar resolution. One may easily end up classifying the 
combined SED of a prestellar core and a nearby, more 
evolved YSO as a Class protostar (see, e.g., CB243, 
Fig.|22]D. 
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5. Description of individual globules 

In this section, we describe in somewhat more detail 
the individual globules, their morphology and physical 
properties, and discuss the evolutionary stages of the 
embedded sources. Distance references are not given 
in the text, but are compiled in Table[2] Globule sizes 
and morphology are summarized in Table|6]and prop- 
erties and evolutionary stages of the embedded sources 
are summarized in Table|7] 

CB4 is a small, roundish, and relatively isolated 
dark globule at a somewhat uncertain distance. Al- 
though |Wood et al.| (Q994) assumed a distance of 
200 pc to the Cas A dark clouds in this region of the 



sky, we follow the analysis of Dickman & Clemens 



1983[ l who derived a distance of 600 pc. This latter 
value is also close to the distance of the "-12kms 



clouds (800 pc; Dame et"aL]|1987| l with which |Laun- 
hardt & Henning associated the globules CB 4 



and CB 6. CB 4 is associated with an IRAS point 
source, located at its southern rim, which is detected 
at lOOyum only. There is no further evidence of an 
embedded source, including the Spitzer maps up to 
160jum. Deep IRAS aperture photometry and molecu- 
lar line studies showed that the globule is cold (~ 7 K) 



and has a relatively low mean density (Clemens et al. 



T99T)|Kane et al.|1994[|Turner et al.|1997| l. Despite its 
compact and dark appearance it thus rather resembles 
a cirrus cloud. CCD polarimetry of background stars 
revealed a faint and well-ordered magnetic field with 
a position angle (RA.) of 71° ( jSen et al.|2005| l, which 
is apparently coupled to the surrounding Galactic field 
(Kane et al. 1995| l. Our non-detection of compact ther- 
mal dust emission confirms the assumption that this 
globule does not seem to form stars. 

CB 6 (LBN 613; Fig.[T]) is a small, cometary-shaped 
globule with a long diffuse tail, visible as scattered 
light cloudshine. It is located only a few degrees from 
CB 4 and has the same radial velocity. We therefore 
assign CB 6 the same distance (600 pc) as CB4. At its 
centre, the globule harbours an embedded YSO with 
circumstellar dust emission and a conical reflection 
nebula (RN0 3; Cohen 1980. |Eiroa et al.|1994| l, likely 
representing scattered light in an outflow cone from the 
embedded YSO. Indeed, broad line wings in CO (2-1) 
^Clemens jni][T99TJ and CS(2-1) (Launhardt eTaL] 
[1998| ) indicate the presence of a molecular outflow 
from the YSO, but no outflow maps exist yet. Scap- 
pini & Codella ( 1996| l searched for, but did not detect 



controversial: while [Clemens & Barvainis ( |1988| l de- 
scribed it as Class III YSO, ' Launhardt et al.| ( fT997 l) de- 
rived a Lboi / Lsintn ratio indicative of a Class protostar 
(Andre et al. 1993| l. Our (sub)mm continuum maps 
CFig[T]l show that the thermal dust emission is domi- 
nated by a compact source centered at the origin of the 
NIR nebula, indicating that the bulk of the mm dust 
emission does not originate from an optically thick en- 
velope, but rather from a disk. The disk and the more 
extended envelope probably obscure the southern out- 
flow lobe which is presumably pointed away from us. 
This morphology, combined with the bolometric tem- 
perature of Tboi ~ 180K let us conclude that this is a 
Class I YSO. 

CB6N is a small globule located xlQ' north of 
CB 6. It has no IR source associated and we did not 
detect it at mm wavelengths, indicating that it does not 
contain a star-forming core. 

CB17 (L 1389; Fig.|2]) is a smafl and slightly 
cometary-shaped globule, located near Perseus and 



associated with the Lindblad ring ( |Lindblad et al. 
|1973[ ). It has a long diffuse tail, pointing north-east 
and visible as cloudshine reflection. The distance of 



CB 17 was derived by [Launhardt & Henning ( 1997 1 



via association in projected space and radial veloc- 
ity with other Lindblad Ring clouds. In the direc- 
tion of CB 17, the Lindblad Ring structures have a 
mean distance from the Sun of ^ 300 pc ( [Dame et aL] 
1987] ). HD 25347, a bright GSIH star at a distance 
of 210 + 40pc ( van Leeuwen||20(j7 1 is located about 
-11' (0.65 pc at 200 pc) south of CB 17 and could be 
responsible for the cometary shape and diffuse cloud- 
shine from the rim and tail of CB 17 if it is located 
at the same distance as the globule. Combining the 
possible associations of CB 17 with both the Lindblad 
Ring and HD 25347, and the distances and uncertain- 
ties involved, we adopt a distance of 250+50 pc for 
CB 17. 

The roundish cloud core of CB 17 is associated 
with a dense submm core and a faint and cold IRAS 
point source (IRAS 04005H-5647) that is detected only 
at 60 and lOOjum. CB 17 has been studied extensively 
by various authors using different mm molecular line 
transitions. The core was found to have a mean kinetic 
gas temperature of Ti^i^ ~ lOK and a virial mass of 



3 Mo (e.g., Lemme et al. 1996 CaseUi et al. [2002 



scaled to D = 250 pc). Pavlyuchenkov et al. (|2006 1 



NH3 and C2S. The evolutionary stage of the YSO is 



compared a number of different molecular line maps 
of CB 17 with chemodynamical models. They derive 
a chemical age of x2 Myr, as well as the relative rota- 
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tional, thermal and turbulent energies, concluding that 
the core will most likely fragment. 

Indeed, our deep 1.3 mm continuum map (Fig.|2]) 
shows a double core with 14" separation and a com- 
mon envelope, located at the south-western rim of the 
globule head. The low-SNR 850 A^m SCUBA map 
can, however, not confirm the double-core nature, 
and the source remains undetected at 450fj.m. The 
Spitzer IRAC maps show a 8 yum shadow that agrees 



in position and shape with the double core (Fig. 28 1 



No emission is detected at wavelengths shorter than 
160jum, indicating these cores have no, or only very 
low-luminosity internal heating sources. Since the 
double core is resolved only in the 1.3 mm map, we 
compile the SED only for the combined core, which 
we call SMM. However, a faint point-hke NIR/MIR 
source (IRS), that is clearly not directly detected in 
the (sub)mm continuum, is located ^IQ" northwest 
of SMM and just outside the 8 |/m shadow, indi- 
cates the presence of a third, more evolved source 



in CB 17 (Fig. 28 1. This source dominates the emis- 



sion at wavelengths shorter 160yum, while SMM dom- 
inates at longer wavelengths. Although the MIPS 3 
map at 160jum does not resolve the two positions, 
a shght shift of the peak position from IRS towards 
SMM indicates detectable emission from SMM at this 
wavelength. For compiling the SEDs, we assigned 
15% of the total 160 /im flux to SMM, but our re- 
sults depend only weakly on the adopted flux splitting. 
The morphology and SED of SMM (Tboi ~ 15 K, 
^smm / i'bol > 30%) suggest this is a prestellar (dou- 



ble) core (cf. Caselh et al. 2002 Pavlyuchenkov 



|et al.||2006) l. Source IRS, which is located at a pro- 
jected separation of ~ 2000 AU from SMM2 within 
the same globule core, is tentatively classified as a 
low-luminosity Class I YSO protostar (Tboi > 55 K, 
ismm / i'bol < 2.6%), although an outflow remains to 
be detected (Chen et al., in prep.). A faint and diffuse 
NIR nebula at the position of SMM2 is speculated to 
be related to IRS rather than to a source embedded in 
SMM2, although the nature of a faint compact NIR 
knot at the peak position of SMM2 remains unknown. 
We are currently analyzing interferometric molecular 
line maps, attempting to better constrain evolutionary 
stages and the kinematic relation of the three sources 
(Schmalzl et al., in prep.). CB 17 is thus one of those 
globules that contain two or more sources of differ- 
ent evolutionary stage with a few thousand AU (see 
Sect.|4!2]and Table|6]). 

CB17SW (L1388) is a small elliptical globule 



of similar size as CB 17, located a; 12' south-west of 
CB 17. It has no IR source associated and we did not 
detect it at mm wavelengths, indicating that it does not 
contain a star-forming core. 

CB 26 (L 1439: Fig.|3]l is a small, cometary-shaped, 
double-core Bok globule located ~ 10° north of the 
Taurus-Auriga dark cloud, at a distance of ~ 140 pc. 
An IRAS point source and dense core with signa- 
tures of star formation is located at the south-west tip 



of the western globule core (Launhardt & Henning 



19971. While the eastern globule core (at 4:56:20, 



52:01:30, B1950, in Fig.[3]) is visible as 8yum shadow 
in the Spitzer IRAC maps, the IR source is located 
»4' further west. The eastern globule core was not 
searched for mm emission, but the Spitzer maps also 
do not indicate the presence of an embedded source. 
OVRO observations at the IRAS position of the mm 
dust continuum emission and of the '^CO(l-O) line 
have revealed a nearly edge-on circumstellar disk of 



radius 200 AU with Keplerian rotation (Launhardt & 



I Sargent 2001 1, surrounding a very young (obscured) 
low-mass TTauri star. It is associated with a small 
bipolar near-infrared (NIR) nebula bisected by a dark 
extinction lane at the position and orientation of the 
edge-on disk (Stecklum et al. 2004 1. The star and 
disk are surrounded by an optically thin asymmetric 
envelope with a well-ordered magnetic field directed 
along RA. ~ 25° ( [Henning et al.|200T] l. Furthermore, 
an Herbig-Haro object (HH494) was identified by Ho- 
and S[II] narrow-band imaging, 6. 15 arcmin northwest 
of CB 26 at RA. = 145° ( |Stecklum et ai:]|2004| . The 
HH object is thus perfectly aligned with the symmetry 
axis of the disk and the bipolar nebula. Recently, |Laun-| 
|hardt et aL] P009| l reported the discovery of a small 
coUimated molecular outflow along the same direc- 
tion, that appears to rotate about its polar axis. It is 
speculated that this rotation is related to the possible 
binarity of the central star We are currently analyzing 
higher resolution interferometric CO maps, attempting 
to reveal the mechanism that drives this outflow rota- 
tion (Launhardt et al., in prep.). Due to the edge-on 
configuration, the bolometric luminosity listed in Ta- 
ble|7]is only a lower limit, while the Lsmm / ^^boi ratio 
is an upper limit. Although the IR spectral index in- 
dicates a Class II source, the other evolutionary tracers 



and the estimated age of 1 Myr (Launhardt et al. 2009 1 
suggest it is a Class I YSO. 

CB 34 (Fig.[4]| has the appearance of a Bok glob- 
ule, but, with a distance of ~ 1.5 kpc, it is rather atyp- 
ical and more massive than other Bok globules. It 
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has three dense cores and is associated with numer- 
ous young stars that seem to have formed from this 
cloud (jLaunhardt & Henning|1997[|Alves et al.[1997 



(2006 1 detected a water maser. Polarimetry of back- 



IHuard et al.|2000| ). A water maser ( |G6mez et al.[2006| 
as well as several collimated outflows testify ongoing 
star formation activity. A bipolar CO outflow from 
cool gas ( |Yun & Clemens] |1994b| ), shocked optical 
knots of atomic emission lines from radiative shocks 



(Herbig-Haro objects), as well as H2 infrared jets ( Yun 



& Clemens 1994a Khanzadyan et al.|2002i ) have been 



reported. The chemi cal age of the globule was esti 
mated to be > 10^ yr ( |Codella & Scappini|l998| . The 
presence of a pre-main-sequence star, CB 34FU, with 
an age of ~ 10^ yr supports this young age ^Alves et al. 



1997] ). [Codella & Scappini| ( [2003| studied the molecu 



lar gas of this globule and showed that the current star- 
formation activity is concentrated in the three main 
clumps which have sizes of ~ 0.25 pc. Because of this 
complexity, we do not further investigate the individ- 
ual sources in this globule, but only present our data. 

CB 39 is a small and not very opaque globule, lo- 
cated in direction of the Gemini OBI cloud. It is 
associated with the well-known Herbig Ae/Be star 
HD 250550, the distance of which was derived with 
700 pc. No mm dust continuum emission has been de- 
tected from this star, indicating that its circumstellar 
disk, if still present, cannot be very massive. Yun &| 
|Clemens| ( | 1 994b| l found a small, bipolar, low- velocity 
(~ 3kms"^) CO outflow with not well-separated 
lobes. CB 39 has been observed in different molec- 



ular lines. It was detected in SO by Codella & Mud- 
1997[l but not detected in CS by |Launhardt et al. 



ers 



( 1998| l. Sen et al.] (2005) measured a magnetic field 
P.A. of 150°, based on the polarization of a number of 
background stars along the line of sight of the globule. 

CB54 (LBN 1042; Fig.|5]l is a relatively compact 
molecular cloud with the appearance of a Bok globule, 
but seems to be associated with the Vela OBI cloud 
complex at a distance of 1.1 kpc. It contains a clus- 
ter of embedded young stars, visible on NIR images 
dYun & Clemens] 1994a||Yun| 1996 '), as well as a cold 
IRAS point source (07020-1618) at the position of the 
mm peak and the embedded star cluster. CO obser- 
vations revealed the presence of a collimated bipo- 
lar molecular outflow with well-separated lobes and 



centered on the IRAS point source (Yun & Clemens 
|1994b). This source is also ass ociated with cm contin- 
uum emission (VLA AY0073; |Yun et al.|l"996j . C'^^O 
observations suggests the presence of multiple cores 
withinfall signatures ( Wang et al.|1995| l. Gomez et al. 



ground starlight revealed a magnetic field with an av- 
erage P.A. of 1 16° ( |Sen et al.|[2005] l. Because of this 
complexity, we do not further investigate the individ- 
ual sources in this globule, but only present our data. 

CB58 (Fig.|6| is a relatively large, slightly elon- 
gated and cometary-shaped globule at a distance of 
1.5 kpc. It is associated with IRAS point source 07161- 
2336 and the NIR images indicate the presence of 
several embedded YSOs. [Afonso et al.| ( |1998| l con- 
cluded from HCN observations that the core repre- 
sents a Class 1 evolutionary stage. Our (low-quality) 
mm dust continuum map suggests the presence of two 
dense cores in the head of the globule, where also the 
IRAS source and a number of reddened stars with as- 
sociated nebulosities are located. However, none of 
the IR sources can be directly associated with one of 
the two dust clumps. Because of this complexity, we 
do not further investigate the individual sources in this 
globule, but only present our data. 

BHR12 (CG30, DC 253.3-1.6; Fig.|7]) is a small, 
opaque, and very pronounced cometary globule with 
a bright rim, located in the Gum nebula region. The 
distance towards the Gum nebula clouds is somewhat 



uncertain, with estimates ranging from 200 pc ( Knude 



etal. 1999 ) to450pc (Reipuith 1983). For consistency 



with earlier papers (e.g., ,Henning & Launhardt|1998 1, 
we use here 400 pc. The globule has an angular diam- 
eter of ~ 2.3 ' and its tail (which merges with the tail 
of CG 31) has a length of 25 ' ( |Reipurth|1983] l. At its 
head side, CG 30 harbors two protostellar cores with a 
projected separation of ~20" (~8000AU), embedded 
in a common envelope ( |Launhardt et al.|[2000j ). The 
northern source (SMM 1, CG30N) is associated with 
the IRAS point source 08076-3556 and a NIR source. 



which drives the Herbig-Haro flow HH 120 (see Ho 
dapp & Ladd 1995 and references therein). The 
southern core (SMM 2, CG 30S;' Launhardt et al.|200l] ) 
is the origin of a protostellar jet with position angle 
44° ( [Hodapp & Laddl[T995) ), but no NIR counterpart 
is seen at this position. Both cores are detected in all 
IRAC bands ( |Chen et al.||2008b| l. The IRAC 4.5 fim 
image also shows bipolar jets, which have been traced 
in CO(2-l) using the SMA ( jChen et al.|2(j08^ . Fig.|7] 
shows that the two cores are well resolved at 450 and 
850 jum, though they are surrounded by a common en- 
velope. The cores are unresolved in our 1.3 mm map. 
The SED is well sampled for both cores, and shows 
that SMM 2 is colder or more heavily extincted than 
SMM 1. We classify SMM 1 as a Class I YSO, based 
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on its bolometric temperature (see Table[7]l. The low 
bolometric temperature of 50 K suggests that SMM 2 
is a Class protostar and thus younger than SMM 2. 
BHR 12 is thus one of those globules that contain two 
or more sources of different evolutionary stage within 
a few thousand AU (see Sect. |4.2| and Table[6]|. 

BHR36 (DC 267.4-7.5; Fig.jS) is a cometary- 
shaped globule in the Vela sheet at a distance of 400 pc. 
A dense core, mapped at 1.3 mm continuum, is located 
at the head side of the globule. The mm emission peaks 
at the position of the IR source detected both by IRAS 
and Spitzer and at the center of a bipolar jet seen on 
NIR images. Two Herbig-Haro objects, HH46/47, are 
associated with a collimated bipolar outflow centered 
at this position (Dopita et al., 1982j. From their ob- 
servations of '^CO(2-l), |Lohr et al.] ( |2007| l derive a 
total mass of 255 + 60 M© and a size of ~ 1.3pc for 
the entire globule. In their ammonia survey, |Bourke| 
|et al.| ( |1995b| l obtain a mass of ~ 11 Mq within the 
central 0.3 pc for the dense core. They also derive a 
kinetic temperature of 7"^ = 18K in the centre. The 
chemical age is of the globule was estimated to be 
~ 0.1 Myr ( [Cecchi-Pestellini et al.|[200T] l. BHR 36 is 
too far south for the JCMT to observe. Like for the 
other southern globules, we are therefore limited to 
the 1.3 mm maps. The fit to the SED of this source 
indicates Tboi of a!l50K (see Table[7]i. We therefore 
classify the embedded YSO as a Class 1 source. 

BHR 55 (DC275H-1.9; Fig.|9} is an average-sized, 
low-extinction globule, associated with the annulus of 



cometary globules in the Vela-Gum complex ( Bourke 



et al. 1995a I, though it has not a very pronounced 
cometary appearance. [Bourke et al.| ( |1995b) l estimate 
the distance to the globule to be 300 pc by looking 
at the reddening of background stars. From ammo- 
nia observations, Bourke et al.| ( |1995b[ ) derive kinetic 
and rotational temperatures of 13K and 12 K respec- 
tively, a central volume density of 8.1 x 10^ cm"-', and 
a core mass of 3.7 Mq. The globule is too far south 
to observe with JCMT, therefore only the 1 .3 mm data 
constrain the long-wavelength end of the SED. The 
very compact (unresolved) mm source is associated 
with an IRAS point source (09449-5052), but has no 
near-infrared counterpart. Due to the lack of data at 
submm and MIR wavelengths, the SED is not well- 
constrained. However, the low Tboi of ai40 K suggests 
this is a Class protostar 



BHR 71 (DC 297.7-2.8; Fig.[TO]l is a very opaque, 
elongated globule located near the southern Coalsack 
at an estimated distance of ~ 200 pc. At its center, the 



globule contains a dense protostellar core, which emits 
strong mm dust continuum emission and is detected by 
various instruments down to MIR wavelengths. Ob- 
servations with higher angular resolution revealed that 
this is actually a double core ([Myers & Mardones 



[T998| ) with two embedded IR sources, IRS 1 and IRS 2 
separated by 17" (3400 AU - [Bourke [[200T] l. Out- 
flows are seen in the H2V- 1-2S(1) line emanat- 
ing from both sources. A highly collimated molecular 
outflow has been detected ( [Bourke et al.|1997] l which 
is associated with two Herbig-Haro objects (HH320 
and HH321; [Corporon & Reipurth[ [1997) 1. jBourkej 
pOOlj l and [Parise et alT] ( |2006| l resolved it into two 
outflows, mapping CO and methanol, and found that 
the IRS 2 outflow is significantly hotter tha n the IRS 1 
outflow. From '^CO(2-l) measurements, Lohr et al. 
( |2007| l derived a mass of 12 + 4 M© for the entire 
cloud. From their ammonia measurements, |Bourke| 
jet al.| ( [T995b| l derived a mass of 2.2 Mg for the dense 
core and 'Bou rke et al.| ( |1997| l derive 2.4 M© from the 
1.3 mm dust continuum map. The 1.3 mm contin- 
uum image does not resolve the two sources, but the 
low Tboi of the combined SED indicates that at least 
IRS 1 (the brighter source at the mm peak position) is 
a Class protostar. This is in agreement with earlier 
work ( jBourkeet al.|1997| l. 

^ cometary- 



BHR86 (DC 303.8-14.2; Fig. 11 



shaped globule located in the Chameleon dark cloud 
complex at a distance of 150 to 180 pc. |L6hr et al. 



(2007 1 carried out '^CO(2-l) observations and mea 
sured a size of 15 ' x 15 ' and a mass of 19 + 4 Mq. 
From ISO observations, jLehtinen et aL] ( [2005| l derived 
central dust temperature of 14 + 1 K and a mass of 
2.7 Mq for the dense core and concluded that the em- 
bedded source is a Class O/I transition object. They 
also showed that the main heating source of the dense 
core is not the embedded source, but the interstellar ra- 
diation field. From CS(2-1), HCO+(1-0) and HCN(1- 
0) fine profiles [Lehtinen| ( [T997] l could classify BHR 86 
as a core with simultaneous infall of dense gas in the 
centre and a non-collapsing envelope. They also de- 
tected a bipolar molecular outflow in '^CO(l-O). NIR 
observations in JHKs showed an exponential extinc- 
tion profile with p - 2.29 ± 0.08 (Kainu lainen et al] 
[2007 ). From the low bolometric temperature of *60 K, 
the high Lsmm / i'boi ratio of s!6%, and the absence of a 
NIR source we conclude this this is a Class protostar, 
which is in agreement with the ISO observations. 
CB 68 (L 146; Fig.[T2| is a small, nearby, slightly 



cometary- shaped opaque Bok globule located in the 
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outskirts of the pOph dark cloud complex at a dis- 
tance of ~ 160 pc. The dense core of the globule is 
associated with a cold IRAS point source (16544- 
1604) and exhibits strong, extended, but centrally 
peaked sub(mm) dust continuum emissio n ([Launhardt 



& Henning||199 7t|Lau nhardt et al.|[T998HVallee et al. 



2000 1. No NIR source is detected at 2.2 jum. The 



central source, which was classified as a Class proto- 
star ( Launhardt & Henning| 19971 1, drives a weak, but 
strongly colhmated bipolar molecular outflow at RA. 
142° fWu et al.|[T996l [M ardone s et al.|[T997l |VaUee| 
|et al. 2000). Vallee et al. (2000) also performed submm 
continuum polarization measurements at SSO/zm and 
found a large linear polarization of 10%, indicating 
that Bok globules can have well ordered magnetic 
fields, and suitable conditions for grain alignment. 
The SED, which is well-sampled between the 1.3 mm 
IRAM data and the NIR IRAC data, indicates a Tboi 
of ~ 50 K (Table[7]i, confirming the classification as a 
Class protostar (Fig.[26|. 

BHR137 (DC 344.6-04.3; Fig.[T3]l is a small, 
opaque, and elongated globule with a diffuse tail. 
The distance is somewhat uncertain. IVilas-Boas et al.l 



(2000j) associate the core with the nearby Scorpius re 
gion at a distance of 145 pc, whereas ( [Bourke et al. 



1995b| l derive a distance of 700 pc, based on measur- 



ing the reddening of background stars. We adopt here 
the latter distance estimate. At the northern rim of 
the globule, outside the dark core, there are several 
red stars with associated nebulosities. An IRAS point 
source (17181-4405) appears to be associated with 
the southern-most of these stars. For the dense core, 
which is located south of the nebulous stars and elon- 
gated north-south, |Bourke et al.| ( |1995b| l estimate a 
lower mass limit 8.7 Mg from NH3. We derive a total 
mass of 26 Mq from the 1.3 mm dust continuum map. 
BHR 137 was not observed by Spitzer and our obser- 
vations do not permit to compile a reliable SED for 
this source. Therefore, we cannot evaluate the evolu- 
tionary stage of this core, although the high mass and 
lack of NIR sources at the position of the mm peak 
suggest a rather early evolutionary stage. Due to the 
large distance and low angular resolution of our SEST 
map, we can also not evaluate the possible multiplicity 
of this massive elongated dust core. 

CB130 (L507; Fig.[T4| is a smafl, slightly elon- 



gated globule located in the Aquila rift at a distance of 
200 pc. There is no IRAS point source associated with 
this globule. [Launhardt & Henning| ( |1997] l detected 
only weak mm-continuum emission in its centre and 



|Yun et aL] ( |1999| l concluded from the weakness of the 
HCN(l-O) emission that there is no embedded YSO. 
Since there is also no evidence of outflows, the source 
was originally classified as a pre-stellar core by |Laun- 



hardt & Henning ( 1997 1. However, our (sub)mm maps 



show that CB 130 is well detected in all bands (see 
Fig. [14]). The core is resolved into two separate com- 
pact sources, separated by 30" (~6000AU) and sur- 
rounded by a common envelope. The brighter source, 
SMM 1, is associated with a faint NIR source, which 
is also detected by Spitzer up to 24 pm (there are no 
useable data at longer wavelengths). A very red star, 
which we caU IRS, is located 14" (a!3000AU) east 
of SMM 1 (Figs.[T4]3 and|28]l. At 24 /im, this star is 



still 1.3 times brighter than the IR source associated 
with SMM 1 . However, as none of the (sub)mm maps 
shows a local maximum at this position, source IRS 
does not seem to have detectable mm dust continuum 
emission associated. We conclude that the spatial co- 
incidence between IRS and the dust core is a projec- 
tion effect and the source is located outside the dust 
core. We could compile a SED for SMMl, sepa- 
rating the fluxes from IRS at short wavelengths and 
from SMM 2 at (sub)mm wavelengths. The low Tboi 
of »55 K and high Lsmm / ^'boi ratio of a; 10% suggest 
SMM 1 is a Class protostar, although the NIR/MIR 



SED is rather indicative of a Class 1 YSO (see Fig. 25 1. 
The IR SED of IRS (aiR = -0.7) suggests this is 
a Class 11 YSO, but we don't have strong constraints 
on its FIR to (sub)mm SED. The (sub)mm SED of 
SMM 2 has a somewhat steeper slope than SMM 1 and 
there is no IR counterpart detected, indicating that this 
core does not have an embedded heating source. We 
therefore assume that this is a very low-mass prestellar 
core. CB 130 is thus one of those globules that contain 
two or more sources with different evolutionary stage 
within a few thousand AU (see Sect. |4.2| and Table[6|. 



CB 188 (L 673-1; Fig.[T5]) is a smafl, roundish, and 
opaque globule with a detached diffuse tail. It is lo- 
cated towards Aquila at a distance of ~ 300 pc and 
is associated with the Lindblad ring. The dense core 
of the globule is associated with a cold IRAS point 
source (19179H-1129) and exhibits strong, extended, 
but centrally peaked sub(mm) dust continuum emis- 
sion. At the peak of the mm emission and the IRAS 
position, there is a very red star, indicating that the 
bulk of the compact mm dust emission does not arise 
from an optically thick envelope, but rather from a cir- 
cumstellar disk. The source drives a bipolar molec- 
ular outflow at RA.= 75° with a mass loss rate of 
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~ 6.6 X 10"^ Mo/yr ( |Yun & Clemens||l994b| ). The 
red and blue lobe have a strong overlap, indicating that 
the source is seen almost pole-on. This is consistent 
with seeing a strong NIR source at the peak of the 
mm continuum emission. The source shows only weak 
N2H^(l-0) emission which can be caused by destruc- 
tion due to CO, which is released from grain surfaces 
in outflows ( |Chen et al.|[2007l l. Our fit to the SED, 
which is well-sampled between 1 .25 fim and 3 mm, in- 
dicates Tboi -300 K. We conclude that CB 188 SMM is 
a Class I YSO, exhibiting a relatively massive circum- 
stellar disk, surrounded by an extended envelope. We 
thus confirm earlier classifications (e.g., |Afonso et al] 
[19981 ). 



CB 199 (B 335, L663; Fig.[T6]i is a roundish, very 
opaque, relatively isolated, and well-studied prototyp- 
ical Bok globule, located towards Aquila (see |Stutz| 
|et al.t2008, and references therein). It's distance was 



recently redetermined to be only 90 - 120 pc ( Olofsson 



& Olofsson|2009[ l. The globule contains a single com- 



pact core at its center, that is detected at far-infrared 
wavelengths and is very bright in the submillimetre 
range (Keene et al.p983| [Chandler et al.|1990l l. [aiou[ 
[et al. (1993 ) have found evidence for mass infall onto 
the central source, which was confirmed by observa- 
tions in a CO rare isotope (Chandler & Sargent' 1993 1 
and by CSS observations (Velusamy et al. 1995j ). The 
source drives a well-coUimated bipolar molecular out- 
flow ( Frerking & Langer|1982 1. Furthermore, a num- 
ber of HH objects have been detected both at optical 



wavelengflis ( |Vrba et al.[[1986[ Reipurdi et aLl[1992[ ) 
and in the 1-0 S(l) fine of H2 (Hodapp[ [1998] l. This 
combination of characteristics that makes the central 
embedded source in B 335 one of the clearest exam- 



ples of a Class protostar. Gomez et al. ( 2006 1 found 
a water maser with a high velocity shift (~ 30km/s) 
compared to the cloud velocity. This could be evidence 
for a protoplanetary or young planetary disk. The low 
Tboi of 37 K that is derived from our fit to the SED cor- 
roborates the earlier classification. 



CB205 (L810; Fig. 17 1 has the appearance of a 
Bok globule, but, with a distance of 2 kpc, its appar- 
ent size of 8 ' corresponds to a diameter of 4.8 pc. It 
is thus much larger and more massive than other Bok 
globules. SCUBA observations revealed the presence 
of a large submm core with two condensations, sep- 
arated by s; 45 " ([Huard et al. 2000'). The first and 



more massive core (CB205-SMM1) coincides with the 
IRAS source 19433H-2743 and has a mass of 10 Mq. 
The total of mass of the dense core is estimated to be 



120 Mq. The cloud is known to harbour many optical 
and NIR sources. ' Neckel & Staude| ( [r990| l showed that 
at least 12 stars are associated with CB 205 and are lo- 
cated in the eastern part of the cloud. A very bright 
and red source (L 8 10 IRS) with a bolometric luminos- 
ity of 890 Lq is the illumination source for a NIR neb- 
ula ( |Yun et al.|[T993) [Massi et al.|[2004] l. This nebula 
is hourglass-shaped and has an extent in NS direction 
of 20". It has the same orientation as the molecular 
outflow, which shows a strong overlap of the blue and 
red lobes. The peak of '^CO intensity map coincides 
with the IRAS source ( |Yun & Clemens|[T994bl l. In 
the western part, [Codella et al. ( 2006| l detected three 
Class sources with enhanced high velocity outflow 
activity, suggesting that the star formation propagates 
through CB 205 as the NIR cluster represents a later 
evolutionary stage. Neckel et al. ( 1985 1 also found 



water maser emission of 7.5 Jy towards L810, which 
could not be confirmed in later observations ( [Gomez | 
et al.|2006 1. Although the entire source was classified 



as Class I source by Afonso et al. ( 1998), we find that 
this complexity does not allow us to construct mean- 
ingful SEDs or to reliably interpret a possible inte- 
grated SED. Therefore, we do not further investigate 
the individual sources in this globule or discuss its evo- 
lutionary stage, but only present our data. 



CB 224 (L 1 100; Fig.[T8]) is a smafl, roundish glob- 
ule with a large visible cloudshine halo (like most 
globules). It is located in the northern Cygnus re- 
gion at d ~ 400 pc. The core of the globule con- 
tains two mm sources, separated by =s35", both de- 
tected at all three (sub)mm bands. The northeast- 
ern source (SMM 2) is associated with the cold IRAS 
source 20355H-6343 and with a red star which is also 
seen at optical wavelengths. The SED of this source 
indicates a Tboi of 340 K and we classify this source 
as a Class I YSO. The southwestern source (SMM 1), 
which is the brighter mm source, is associated with a 
very faint and diffuse NIR nebula. This source drives 
a colhmated "CO bipolar outflow at RA. = -120° 
( [Chen et al.[[2007] ). HCO+(1-0) observations clearly 
show evidence of infall towards the core center ( |De[ 
[Vries et al"l [2002). CB224 was not observed with 
Spitzer. The SED is therefore not well-constrained at 
MIR/FIR wavelengths and we rely on aperture pho- 



tometry in IRAS maps [Gezari et al.| ( [1999| l. The SED 
of this source indicates a Tboi of 56 K and we classify 
it as a Class protostar. 

CB230 (L1177; Fig.[T9|) is a smafl, slightly 
cometary-shaped, bright-rimmed Bok globule lo- 
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cated in the Cepheus Flare region at a distance of 
~ 400 ± 100 pc. The dark globule has a low-extinction 
tail, which is clearly visible on optical images by its 
cloudshine, that connects to the bright nebula Sh2- 
136. At its southern rim, the globule contains a dense 
core that is associated with two NIR reflection nebu- 
lae, separated by ~ 10" ( | Yun|| 1 996| |Launhardt| 1 996 



2001 ). The brighter, western nebula (IRS 1), has a con- 



ical shape with the origin at the mm peak position. The 
fainter eastern nebula (IRS 2) is also located inside the 
mm core, but seems to exhibit much less dust emission 
than IRS 1. An IRAS point source (21 169H-6804) is lo- 
cated between the two NIR nebula, indicating that both 
sources contribute to the IRAS PSC fluxes. CB230 
was classified by Froebrich ( 2005 |l as a Class O/I tran- 
sition object. CB230 was also observed in several 
molecular lines and was found to exhibit signatures 
of mass infall. Its virial mass was determined to be 



12 Mo (CS) and 8 Mq (C'^O), respectively ( jLaun- 



Ihardt et al.|T998)|Launhardt|2001 ). The dense core is 
associated with a large-scale coUimated CO outflow 
with 6.6 X 10"^ Mp/ yr at RA. - 7° and a dyn amical 
age of ~ 2 X lO'^yr ( |Yun _& Clemens "l994b|). This 
outflow was resolved by ILaunhardt (2001|l into two 
misaligned bipolar flows, the stronger one originating 
at IRS 1 and the weaker one at IRS 2. The NIR neb- 
ulae can thus be interpreted as forward-scattered light 
emanating from those outflow cones that are directed 
towards us. Observations of [Fell] and weak H2 emis- 
sion in the jet suggest the presence of fast, dissociative 
J-shocks (Ma ssi et al.||2004| l. Magnetic field strength 
and projected direction in the dense core were derived 
by submm polarimetry to be B = 218 ±50 jiG and RA. 
= -67°, respectively ( Wolf et al.|2003 ). Earlier work 



by |Vallee et ai:] ( |2000 



is consistent with this picture. 
Since the two embedded sources are only resolved at 
wavelengths shorter 25 /urn, we cannot compile com- 
plete separate SEDs. However, the NIR luminosity 
ratio and the submm maps indicate that IRS 1 is the 
more massive and luminous source that dominates the 
SED at all wavelengths. The combined SED and the 
resulting Tboi of ai210K should thus be representative 
for IRS 1 . While the bolometric temperatures and high 
■^'smm / ibol f^tio suggcsts a Class I YSO with a signifi- 
cant envelope, the aiR value of the double-peaked SED 
indicates Class II (see discussion in Sect. |4.1| l. From 
flie similarity of the NIR/MIR SEDs and the spatial 
proximity, we also conclude that IRS 2 is at the same 
evolutionary stage, but less massive, and has formed 
from the same initial core as IRS 1 . 



CB 232 (B 158; Fig.|20| is a small, sHghtly cometary 
shaped globule which is associated with the Sh 2-126 
cloud at a distance of 600 pc. Using SCUBA obser- 
vations, |Huard et al.| ( |1999| l detected a submm double 
core with the two components having different evo- 
lutionary stages and masses of 0.3 M© and 0.1 M© 
respectively. The total mass was derived to be 6 Mq. 
The more evolved object, which we call here IRS and 
which presumably can be associated with a nearby 
NIR star and IRAS source 21352H-4307, was classified 
to be a Class I source ( |Yun & Clemens||l995| l. The 
main submm core (with no NIR counterpart; called 
here SMM) is a candidate Class protostar. A bipolar 
outflow with poor collimation at R A. - -30° ( Yun & 



Clemens|1994b| l cannot clearly be associated with ei- 
ther of the two sources. The virial mass of the entire 
dense core was derived with 42 Mq (CS) and 16 Mq 



(C"*0), respectively ( jLaunhardt etal.|l998| . CB 232 is 
also currently the smallest Bok globule that is known 
to harbour a water maser ( |G6mez et al.||20()6) l. Al- 
though our submm maps, except for the low SNR 
450 jum map, do not resolve the two embedded sources 
(SMM and IRSl), we tried to decompose the emission 
into an extended envelope and two compact gaussian 
sources at the positions of the mm peak and the red 
star. The resulting individual source fluxes at 850 fj.m 



are somewhat higher than estimated by Huard et al. 



( |1999| l and the relative uncertainty is of the order 
50%, but they can be used to constrain the SED of 
source IRS fairly well. The Spitzer maps show that 
even at 70 fim source IRS 1 completely dominates and 
there is no trace of emission detected from the position 
of source SMM, indicating that the dust core has no 
or only a very low-luminosity internal heating source 
(see Fig. 28 for the 24yum image). Source SMM is 



thus more likely prestellar than protostellar The SED 
for IRSl, which may contain some unknown but small 
contribution from SMM around 100/im, indicates a 
Tboi of * 170 K. We classify IRSl as a Class I YSO. A 
second IR source (IRS2), located »40" south of SMM, 
is seen on all Spitzer images and seems to be associ- 
ated with the globule, but we don't compile an SED or 
classify this source. CB 232 is one of those globules 
that contain two or more sources of different evolu- 
tionary stage within a few thousand AU (see Sect. |4.2| 
and Table|6]). 



CB 240 (L 1 192; Fig.|2T]) is a smafl, roundish glob- 
ule with a diffuse tail, located in the Cepheus Flare 
region. Its distance is estimated to be ~ 500 pc, but 
this is rather uncertain. IRAS source 22317H-5816, 
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which has been identified as a Class II-D YSO ( |Yun &| 
|Clemens|1995[ ), is located outside the visible boundary 
of the globule, =s2' to the south-east. In a NIR survey, 



Yun & Clemens ( 1994a| l detected a red object at the 



IRAS position, but no nebulosity or evidence of mul- 
tiplicity. Launhardt & Henning|(|19971l, however, iden- 
tify a cluster of young stars at the IRAS position (see 



Fig. 21 3), and calculate the bolometric luminosity of 
the associated core to be 7.3 Lq. CB 240 has been de- 
tected in a number of molec ular transitions, including 
C^^O, NH3, HCO+ and SO ( Wang et al.| 1995] [Turner 



|1995a|b|cl l. No outflow has been detected from this 
source ( |Yun & Clemens] 1992"| l. CB 240 is detected in 
our dust continuum maps with a low signal-to-noise ra- 
tio at 850 ijm and 1 .3 mm. It is not detected at 450 fim. 
Since the IRAS source and mm core is located outside 
the optical boundaries of the globule and is obviously 
composed of a small cluster of young stars with a dif- 
fuse common dust envelope, we cannot discriminate 
and classify individual sources. Therefore, we do not 
show a SED and do not further discuss this source here. 
CB 243 (L 1246; Fig. 21 1 is an opaque, very elon- 



gated, filamentary-like globule with a long dark tail 
pointing eastward, located at a distance of 700 pc in the 
CepOB3 region. A head-like, very opaque extension 
at its western end contains a dense core that is currently 
forming one or more stars. Extended sub mm dust 
emission is traced from the entire western part of the 
globule, with two embedded compact sources: SMM 1 
in the western head of the globule and SMM 2 towards 
the center of the globule ( [Visser et al.|2001| l. SMM 1 
is associated with the IRAS source 23228H-6320 and 
was classified a Class 0/1 transition source ( |Froebrich 



2005 1 with a bolometric luminosity of 10 Lq ( [Laun- 



hardt & Henningl[T997] l. [Visser et al.[ ( pOOT l) also de- 
tected a molecular outflow from SMM 1 in '^CO(2-l), 
with a dynamical age of 4.6 x 10"* yr and masses of 
0.004 Mq in the red and 0.015 Mq in the blue lobe, re- 
spectively. Our data confirm the general morphology 
of the dust emission from this globule. However, we 
find that the dust core SMMl is clearly offset by ~ 10" 



to the north-west with respect to IRS (see Fig. 21 ->) 



This offset is larger than the pointing uncertainty of 
the mm telescopes and is confirmed independently in 
both the IRAM and SCUBA maps. Unfortunately, the 
Spitzer IRAC maps cover only the bulk of the glob- 
ule containing SMM 2, but not the head with SMM 1. 
They show that SMM 2 has an 8 jum shadow and does 
not contain an embedded heating source. We there- 
fore assume that SMM 2 is a prestellar core. How- 



ever, the MIPS 24 map shows clearly that there is 
a point-like MIR source at the IRAS position close to 
SMM 1, but no emission is detected at the center posi- 
tion of SMM 1, indicating that the dust core has no or 
only a very low-luminosity internal heating source (see 



Fig. 28 for the 24 //m image). Since we cannot separate 
the SEDs of the two sources (SMM 1 and IRS) and it 
is unclear from which of them the molecular outflow 
originates, we cannot clearly classify the sources in 
terms of their evolutionary stage. We can only spec- 
ulate that IRS is a low to intermediate-mass Class I 
YSO that has most likely formed in CB 243, and that 
SMM 1 is most likely a 5 - 10 Mq prestellar core. This 
of course also makes invalid the classification attempts 
by Visser et al. ( 2002) l and [Froebrich ( 2005[ l, who as- 
sumed that all fluxes arise from one source. CB 243 is 
thus one of those globules that contain two or more 
sources of different evolutionary stage within a few 
thousand AU (see Sect. |4.2| and Table[6]l. 



CB 244 (L 1262; Fig. 23 1 is a medium-sized, slightly 



elongated globule located in the Cepheus region and 
most likely related to the Lindblad ring at a distance 
of » 200 pc. The globule contains two (sub)mm cores, 
surrounded by extended envelopes. The brighter (at 
submm wavelengths) and more compact core (SMM 1) 
is associated with a faint NIR reflection nebula and 
a cold IRAS point source (23238H-7401). The SED 
of this source is well-sampled between NIR and mm 
wavelengths and we derive r^,,; ^ 70 K. The source 
drives a collimated bipolar molecular outflow with 
total mass 0.07 Mq ( [Terebey et aL] [1989} [Yun & 
[Clemens] [1994b ]). We confirm earlier classifications 
of CB 244 SMM 1 as a Class protostar, although it 
is probably more evolved than, e.g., B335. The sec- 
ond core, SMM 2, which is located ^90" (18,000 AU) 
north-west of SMM 1, is much less compact and only 
detected at low signal-to-noise in the 850 jum and 
1.3 mm maps. There is no indication of an embed- 
ded source in this core. Rather, it is seen as a dark 



shadow in the 8yum IRAC A map (Tobin et al. 



and in the 24 //m MIPS map (Stutz et al.][2010 



2010 1 



We 



therefore assume that SMM 2 is a prestellar core (cf. 

2009|l. For the region between the two 



Stutz et al. 



submm-cores. Wolf et aL] ( [2003] ) observed an align- 
ment of the polarization vectors parallel to the gas/dust 
bridge connecting the two cores. This could be a trace 
of mass infall along field lines in the early stages of 
star formation. In the outflow region around SMM 1, 
the polarization vectors tend to align with the direction 
of the outflow. 
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CB246 (L 1253; Fig.|24]l is a medium-sized, dark, 
roundish globule with a narrow tail pointing southeast. 
It is most likely associated with the Lindblad ring, the 
distance to which in this direction is ~ 140pc. The 
globule has an extended complex core that is frag- 
mented into at least two elongated sub-cores, but no 
compact embedded sources have been found. Spitzer 
IRAC images show the two sub-cores as strong 8 yum 
shadows. Lemme et al.| (il996 ) and Codella & Scappini 
( [T998| observed CB246 inl:2S(2-l) and NH3. They 
also identified the two sub-cores and derived their total 
mass with » 7 Mq. They conclude that this source as 
inactive, with neither a water maser nor molecular out- 



flow. Crapsi et al. (20051) concluded from the observed 
low N2D+/N2H+ -fraction that CB 246 is still in a very 



early evolutionary phase. Lee et al. ( 2001) 1 classified 
CB 246 as having significant red excess, but no blue 
asymmetry in the CS spectra, indicating that the core 



has not begun to collapse. Recently, Sen et al. (2005 1 
and Ward-Thompson et al. ( 2009| l derived the orienta- 
tion of the magnetic field by measuring the polariza- 
tion of background stars. They find that the magnetic 
field has a mean angular offset from the short axis of 
the globule core of 28 ± 20 deg. Like for other globule 
they have studied, this hints towards magnetic fields 
not governing the star formation in these clouds. Be- 
cause the source is not detected at wavelengths shorter 
than 450 //m, we cannot compile a meaningful SED. 
We conclude that CB 246 contains at least two star- 
less cores, possibly prestellar in nature, thus confirm- 
ing earlier classifications. Sub-core SMM2 is fairly 
elongated and might be fragmented into a chain of sub- 
clumps. 

6. Summary and conclusions 

We have studied the dense cores of 32 Bok glob- 
ules and have obtained deep NIR images and (sub)mm 
dust continuum maps at up to three wavelengths (0.45, 
0.85, and 1.3 mm). With the exception of a small con- 
trol sample, all sources were selected from earlier sur- 
veys which identified them as good candidates for hav- 
ing ongoing star formation, i.e., this is not an unbi- 
ased survey. We also compiled SEDs, taking special 
care of separating the flux contributions from differ- 
ent neighbouring sources at different wavelengths, and 
fitted them to derive various source quantities and evo- 
lutionary stages of the sources. The main results of 
this study are: 

1. We detected (sub)mm dust continuum cores in 



26 out of the 32 globules observed. The ten- 
tative, low SNR, single-dish 1.3 mm continuum 
detections in CB52 and BHR41, published in 
and earlier paper, could not be confirmed. 

2. Eight of the 26 globules with detected (sub)mm 
cores are not studied in further detail or evalu- 
ated in terms of their evolutionary stage because 
they were either at too large distances (> 1 kpc) 
and multiple embedded sources were not re- 
solved, or the (sub)mm maps were of too low 
quality, or we simply did not have enough data 
to draw any reliable conclusions. 

3. In 18 globules with detected (sub)mm cores, we 
derived evolutionary stages and physical param- 
eters of the embedded sources. In total, we iden- 
tified nine starless cores, presumably prestellar, 
nine Class protostars, and eleven Class I YSOs 
in these 18 globules. 

4. We find that the bolometric temperature is the 
most reliable tracer to discriminate between 
Class protostars and Class 1 YSOs and confirm 
the empirical boundary of 70 K. The spread of 
^smm / iboi ratios within the Class and Class 1 
groups is relatively large (2 - 10% within Class 
and 0.8-3.5% within Class 1), with no signifi- 
cant correlation between Tboi and Lsmm / ^boi 
within the groups. However, the three most 
evolved Class I sources, with visible stars and 
compact (sub)mm emission arising presumably 
from circumstellar disks, also have the lowest 
ismm/i'boi ratios (< 1.3%). We take this as 
tentative indication that the Lsmm / i'boi ratio, 
as indicator of envelope dispersal, may better 
trace the evolution within the Class and Class I 
phases. 

5. The mean FWHM core sizes decrease from 
8,000 AU for prestellar cores, to 3,000 AU for 
Class protostars, to <2,000AU for Class I 
YSOs. The latter group also exhibits extended 
remnant envelopes with diameters of order 
20,000 AU. Source-averaged volume densities, 
n-n, increase from 1 x 10^ cm"^ for prestellar 
cores, to 7 X 10^ cm"-' for Class protostars, to 
> 8 X 10^ cm"^ for Class I YSOs. The extended 
envelopes of the Class 1 YSOs have a mean den- 
sity of 7 X lO'^cm-l 

6. At least two thirds (16 out of 24) of the star- 
forming globules studied here show evidence of 
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forming multiple stars on scales between 1,000 
and 50,000 AU, either as multiple star-forming 
cores, wide embedded binaries, or small star 
clusters. The fraction of closer binaries formed 
from unresolved mm cores might be higher, but 
remains unknown from this study. 

7. We find that the large majority of these small 
prototstar and star groups in globules with mul- 
tiple star formation are comprised of sources 
with different evolutionary stages. This includes 
neighbouring mm sources with obviously dif- 
ferent evolutionary stages, prestellar or proto- 
stellar cores with nearby IR sources, presum- 
ably more evolved protostars or Class 1 YSOs, as 
well as NIR star clusters next to large (sub)mm 
cores with the potential to form more stars. In 
only three globules we find coeval pairs, rang- 
ing from multiple pre-stellar cores in CB 246, 
embedded Class protostars in BHR71, to an 
embedded Class 1 YSOs pair in CB 230. This 
widespread non-coevality possibly suggests a 
picture of slow and sequential star formation in 
isolated globules. 

8. These findings also call for special attention 
when compiling SEDs and attempting to derive 
source properties from flux measurements with 
insufficient angular resolution. One may eas- 
ily end up classifying the combined SED of a 
prestellar core and a nearby, more evolved YSO 
as a Class protostar. 

While this paper presents an extensive, though not 
complete inventory of star-forming cores in nearby 
Bok globules, it remains subject of other detailed stud- 
ies to investigate the physical and chemical proper- 
ties, multiphcity and evolutionaly stages of individual 
sources, and to evaluate if and how star formation in 
isolated globules differs from that in larger molecular 
cloud complexes. We expect soon to obtain spatiaUy 
resolved mid to far-infrared Herschel data for many of 
the globule cores presented here and plan to use these 
data to accurately measure the dust temparure profiles 
and derive more reUable density profiles. In particular 
for the prestellar cores, which characterize the initial 
conditions of the protosteUar coUapse, this may lead to 
significant corrections of previous estimates of mass 
and density distributions which were often hampered 
by the lack of temperature measurements. Further- 
more, we are currently working on detailed studies of 
two globule cores which characterize the stage of the 



onset of the protostellar collapse and which may shed 
hght on the non-coeval evolution of different sub-cores 
within one globule (Schmalzl et al., in prep.). Last, but 
not least, this paper shall provide some guidance for 
further follow-up studies of individual globules, e.g., 
with Herschel, SCUBA2, and ALMA. 
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A. Optical depth conversion for JCMT 

For the JCMT data, the optical depths quoted in Table[3]refers to 225 GHz. To convert to rgso and T450, the following 
conversions ( Archibald et al.|2002[ l can be used for dates prior to October 1999: 

T850 = 3.99(T225Cft " 0.004), (Al) 

T450 = 23.5(r225GHz - 0.012). (A2) 
For dates post December 1999, the following conversions apply: 

T850 = 4.02(T225C//^ " 0.001), (A3) 
T450 = 26.2(T225GH^ " 0.014). (A4) 

The change is due to the different filters which were installed in SCUBA in 1999. 
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Table 1 

Typical properties of Bok globules and their dense cores 





Globule 


Dense Core 


Mass 


5 -50 Mo 


0.5 -5 Mo 


Size 


0.2 - 1 pc 


0.02 - 0.05 pc 


Mean density 


10^ 


10^ cm-3 


Gas temperature 


15 K 


lOK 


Line widths 


0.5-2kms-' 


0.4 -0.7 km S-' 



Table 2 

Coordinates and distances of the observed globules 



Name 


Other names'' 


IRAS source 


R.A." 


Dec.'' 


D (Ref.) 




Observation codes'-' 










(1950) 


(1950) 


[pc] 




NIR 


450&850//m 


1.3 mm 


CB4 






00:36:15.0 


+52:35:00 


600 


(11) 


CA02 




196/97 


CB6 


LBN 613 


00465+5028 


00:46:33.8 


+50:28:25 


600 


(1) 


CA93 


J97a 


195 


CB6N 






00:46:50.0 


+50:39:45 


600 


(1) 






196/97 


CB 17 SW 


L1388 




03:59:16.4 


+56:42:27 


250 


(1) 






196 


CB 17 


L1389 


04005+5647 


04:00:31.3 


+56:47:58 


250 


(1) 


CA95 


J98d 


195 


CB 26 


L1439 


04559+5200 


04:55:54.7 


+52:00:15 


140 


(3) 


CA93 


JOOa 


196 


CB34 




05440+2059 


05:44:04.7 


+20:58:54 


1500 


(2) 


CA93 


J97d/98a 


195 


CB39 


(HD250550) 


05591+1630 


05:59:06.0 


+ 16:30:58 


700 


(4) 


CA93 




195 


CB52 




06464-1650 


06:46:25.3 


-16:50:38 


1500 


(2) 


CA93 




195 


CB 54 


LBN 1042 


07020-1618 


07:02:07.2 


-16:18:58 


1100 


(5) 


CA93 


JOOa 


195 


CB58 




07161-2336 


07:16:08.3 


-23:36:19 


1500 


(2) 


CA93 


J97d 


S02 


BHR12 


DC 253.3-1.6 


08076-3556 


08:07:41.4 


-35:56:12 


400 


(2,6) 


LS95 


JOOa 


S02 


BHR36 


DC 267.4-7.5 


08242-5050 


08:24:16.8 


-50:50:49 


400 


(7,6) 


LS95 


JOOa 


S02 


BHR41 


DC 267.7-7.4 


08261-5100 


08:26:11.5 


-51:00:39 


400 


(7,6) 


LS95 




S02 


BHR55 


DC 275.9+1.9 


09449-5052 


09:44:56.9 


-50:52:07 


300 


(7,6) 


LS95 




S02 


BHR58 


DC 289.3-2.8 


10471-6206 


10:47:07.2 


-62:06:28 


250 


(7,6) 


LS95 




S02 


BHR71 


DC 297.7-2.8 


1 1590-6452 


11:59:03.3 


-64:52:17 


200 


(8,6) 


LS95 




S95 


BHR86 


DC 303.8-14.2 


13036-7644 


13:03:38.9 


-76:44:11 


180 


(9,6) 


LS95 




S02 


BHR87 


DC 307.3+2.9 


13224-5928 


13:22:26.0 


-59:28:11 


400 


(7) 


LS95 




S02 


CB 68 


L146 


16544-1604 


16:54:27.6 


-16:04:48 


160 


(2) 


CA94 


J98b,d 


194/95 


BHR137 


DC 344.6-4.3 


17181-4405 


17:18:08.9 


-44:06:02 


700 


(7,6) 


LS95 




S02 


CB 130 


L507 




18:13:40.0 


-02:34:00 


200 


(2) 


CA02 


JOla 


197 


CB 188 




19179+1129 


19:17:53.9 


+ 11:29:54 


300 


(2) 


CA02 


J97b/03c 


195 


CB 199 


B 335, L 663 


19345+0727 


19:34:35.3 


+07:27:24 


100 


(10) 




J07e/98b,d/01a 




CB205 


L810 


19433+2743 


19:43:22.1 


+27:43:39 


2000 


(2) 


CA93 


J97e/03c 


195 


CB 224 


LI 100 


20355+6343 


20:35:30.1 


+63:42:58 


400 


(5) 


CA93 


J03a,b 


194 


CB 230 


L1177 


21169+6804 


21:16:52.7 


+68:05:09 


400 


(5) 


CA93 


J97b/0Ia 


194 


CB 232 


B 158 


21352+4307 


21:35:14.2 


+43:07:04 


600 


(2) 


CA93 


J98b/00b 


195/97 


CB240 


LI 192 


22317+5816 


22:31:44.1 


+58:16:21 


500 


(2) 


CA93 


J97e/00b 


195 


CB243 


L1246 


23228+6320 


23:22:54.2 


+63:19:54 


700 


(2) 


CA93 


J97b.c/98c 


194/95 


CB244 


L1262 


23238+7401 


23:23:40.0 


+74:01:30 


200 


(12) 


CA93 


J97b,c/98b/01a,b 


195 


CB246 


L1253 




23:54:05.0 


+58:18:00 


140 


(2) 


CA93 


J00b,c/01a 


196/97 



" L= |Lynds|(T962) ; LBN= |Lynds](T965l ; B= |Bamard eTaL](T927) ; BHR= |Clemens & Barvaiiiis](T988) 

''Reference position for all Figures; not always identical to IRAS PSC position 
'^Near-IR observing details, see Table [4| (Sub)mm observing details, see Table[3| 



References. — (1) this paper (Sectjs), (2) Launha'rdt"& Henning'(T99Tl, (3) Launhardt & SargentV2001'l, (4)'Testi et al.T19971, (5)'mif| 
p al. 1 20031, (6) Henning & Launhardt ( 1998 i, (7) Bourke et al. I r995bi. (8) Bou rke et al. (1997 l, (9) Whittet et al. i l997 i, (10) 01ofsson'S| 
Pofison]pXi9) , (ll) |Dickman & Clemens||1983( , (12) |Hilton & Lahulla|)1995| ' ' 
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Table 3 

submillimetre and millimetre continuum observing parameters 



Obs. 


Telescope 


Date 


Detector 




HPBW 


TZenith 


Flux calibrators" 


code 








[mm] 


[arc sec] 







J97a 


JCMT 15m 


06/97 


SCUBA 


0.45/0.85 


8.6 / 14.9'' 


0.062 ± 0.007" 


Ur 


J97b 


JCMT 15m 


08/97 


SCUBA 


0.45/0.85 


8.5 / 14.8 


0.09 ± 0.03" 


Ur, Ma 


J97c 


JCMT 15m 


09/97 


SCUBA 


0.45/0.85 


9.5 / 15.1 


0.105 ±0.008" 


Ur 


J97d 


JCMT 15m 


10/97 


SCUBA 


0.45/0.85 


8.6 / 14.9'' 


0.17 ±0.03" 


Ur, 618 


J97e 


JCMT 15m 


12/97 


SCUBA 


0.45/0.85 


8.6 / 14.9'' 


0.042 ± 0.006" 


HL, 10216 


J98a 


JCMT 15m 


02/98 


SCUBA 


0.45/0.85 


8.6 / 14.9'' 


0.052 ± 0.005" 


10216, 16293 


J98b 


JCMT 15m 


04/98 


SCUBA 


0.45/0.85 


8.2/ 14.6 


0.06 ± 0.04" 


Ur, 10216, 2668,618 


J98c 


JCMT 15m 


07/98 


SCUBA 


0.45/0.85 


8.0 / 14.4 


0.032 ± 0.004" 


Ur 


J9Sd 


JCMT 15m 


08/98 


SCUBA 


0.45/0.85 


8.0 / 14.4 


0.096 ± 0.006" 


Ur, 16293 


JOOa 


JCMT 15m 


03/00 


SCUBA 


0.45/0.85 


8.6 / 14.9'' 


0.28 ± 0.05" 


10216,618, 16293 


JOOb 


JCMT 15m 


06/00 


SCUBA 


0.45/0.85 


8.3 / 14.8 


0.107 ± 0.004" 


Ur, 16293 


JOOc 


JCMT 15m 


09/00 


SCUBA 


0.45/0.85 


8.6 / 14.9'' 


0.088 ± 0.008" 


618 


JOla 


JCMT 15m 


09/01 


SCUBA 


0.45/0.85 


8.5 / 14.9 


0.079 ± 0.009" 


Ur, 618 


JOlb 


JCMT 15m 


10/01 


SCUBA 


0.45/0.85 


8.6 / 14.9'' 


0.056 ± 0.005" 


Ma, 2668 


J03a 


JCMT 15m 


02/03 


SCUBA' 


0.45/0.85 


9.3 / 15.8 


0.042 ±0.001" 


Ur 


J03b 


JCMT 15m 


05/03 


SCUBA' 


0.45/0.85 


8.6 / 14.9'' 


0.037 ±0.001" 


Ur 


J03c 


JCMT 15m 


09/03 


SCUBA 


0.45/0.85 


9.3 / 15.8 


0.080 ± 0.005" 


Ur 


194 


IRAM 3()m 


04/94 


MPIfR 7-chan. 


1.3 


10.5 


0.26 ±0.12 


Ur, Ma 


195 


IRAM 30m 


02/95 


MPIfR 7-chan. 


1.3 


10.5 


0.23 ±0.11 


Ur, Ma 


196 


IRAM 30m 


02/96 


MPIfR 19-chan. 


1.3 


10.5 


0.20 ±0.10 


Ur 


197 


IRAM 30m 


03/97 


MPIfR 19-chan. 


1.3 


10.5 


0.22 ±0.12 


Ur 


S95 


SEST 15m 


11/95 


MPIfR 1-chan. 


1.3 


23 


0.19 ±0.06 


Ur 


S02 


SEST 15m 


1 1/02 


SIMBA 


1.3 


24 


0.20 ± 0.05 


Ur, Np 



^Planets and standard calibration sources. Ur = Uranus, Ma = Mars, Np = Neptune, 618 = CRL 618, HL = HL Tau, 
10216 = IRC+ 10216,16293 = IRAS 1 6293-2422, 2668 = CRL 2668 

''Zenith optical depth at 225 GHz, monitored by the Caltech Submillimeter Observatory radiometer. See Sect.[A]for 450- and 
850~//m conversions 

'^SCUBA operated in scan-map mode 

^Planet beam map not available. Quoted beam-sizes are an average 



Table 4 

NIR observing PARAMETERS 



Obs. code Telescope Date Camera Pixel scale Filters" Seeing 



CA93 


Calar Alto 3.5m 


12/93 


MAGIC 


0'/32 


J, H, KS 


0/9 - 2'/ 2 


CA94 


Calar Alto 3.5m 


08/94 


MAGIC 


0'.'32 


KS 


l'/4 


CA95 


Calar Alto 3.5m 


01/95 


MAGIC 


0'/32 


KS 


a'5 - 1" 


CA02 


Calar Alto 3.5m 


10/02 


nPrime 


0-/40 


J, H, KS 


l'/4 


LS95 


La Silla 2.2m 


03/95 


IRAC-2B 


(X'50 


J, H, KP 


1'/ 1 - l'/4 



"Ao/AAmfim: 1: 1.25/0.30; H: 1.65/0.30; KS: 2.16/0.32; KP: 2.10/0.34 
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Table 5 

Results of the millimetre and submillimetre continuum mapping 



Name 


Peak position" 
(BI950) 
[h:m:s, °:':"] 


450 /jm 

[Jy/beain] 


c lolb 
i3 |, 

[Jy] 


850//m- 

(*^b) 

[Jy/beam] 


t- lolb 
[Jy] 


1.3 mm — 

ipeiik ,f~. 

ti^b) 

[Jy/beam] 


Fig 

^lolb 

[Jy] 


CB4 














<0.013'^ (10.5) 






CB 6N 














<0.014'^ (10.5) 






V 1 J I / fcj vv 














<-n CiO'?^ no 51 






CB 39 














<0.027' (10.5) 






CB52 














<0.023' (10.5) 






BHR41 














<0.060'^ (24) 






CB6 


00:46:33.8, 


+50:28:25 


1.0 (10) 


8.0 


0.18 (14.9) 


0.8 


0.058 (10.5) 


0.18 


_ 

1 


CB 17 












1.3 




0.55 


5 


SMMl 


04:00:33.6, 


+56:47:52 


<1.2' (8.0) 




0.17 (15) 




0.039 (10.5) 




5 


SMM2 


04:00:31.5, 


+56:47:58 


<1.2'^ (8.0) 




0.13 (15) 




0.029 (10.5) 




2 


CB 26 


04:55:54.6, 


+52:00:15 


3.2 (8.6) 


6.7 


0.48 (14.9) 


0.6 


0.16 (10.5) 


0.24 




CB 34 






■■■ 


40 




5.5 




2.0 


5 


olVlMl 


05:44:06.4, 


+20:59:38 


1.2 (8.6) 


10 


A /I /I ^ 1 /I A\ 
U.44 (14.y) 


1 .5 


u.uyu (lu.j ) 


0.70 


1 


SMM2 


05:44:02.5, 


+20:59:05 


1.0 (8.6) 


5 


0.29 (14.9) 


0.7 


0.071 (10.5) 


0.40 


1 


SMM3 


jy.y. 


+zU. jy. J 1 


0.5 (8.6) 


2 


0.13 (14.9) 


0.2 


0.045 (10.5) 


0.30 




CB 54 


07:02:06.0, 


-16:18:51 


9.3 (8.6) 


74.0 


1.5 (14.9) 


7.3 


0.34 (10.5) 


2.05 


5 


CB58 


07:16:09.4, 


-23:35:52 






<0.27' (14.9) 




0.12 (10.5) 


0.48 


% 


BHR 12 








27.5 




D.y 


A tlA\ 
U.O/ 1,Z4J 


1.54 


7 


olvilvil 


08:07:40.3 


-35:56:05'' 




j.y 


1.2 (14.9) 


1 "7 




0.49'' 


7 


SMM2 


08:07:40.0 


-35:56:26'' 


3.7 (8.6) 


6.2 


1.0 (14.9) 


1.2 


e 


0.36'^ 


7 


BHR 36 


08:24:15.8, 


-50:50:44 






1.2 (16) 


5.0 


0.2S (24) 


0.85 


S 


BHR 55 


09:44:57.2, 


-50:52:11 










0.19 (24) 


0.42 




BHR 58 


iu.4 / .Ul .u. 


— OZ.UO. jU 










0.13 (24) 


u. iz 




BHR 71 


1 1 .jy.uj.u. 


— 64'5 1 '58 










1.83 (23) 


J.0\) 


10 


BHR 86 


13:03:41.6, 


-76:44:23 










0.26 (24) 


1.33 


TT 


BHR 87 


13:22:24.5, 


-59:27:57 










n 1 1 (0A\ 
u. ii y^^) 


1.00 




CB 68 


16:54:27.2, 


-16:04:44 


2 5 (8 6) 


21 


53 n4 6^ 


4 


U. 1 J 1, lU. J ) 


1.55 


12 


BHR 137 


17:18:08.5, 


-44:06:17 










0.38 (24) 


1.25 


t5 


CB 130 








12 




2.2 




1.20 


t5 


SMMl 


18:13:39.6, 


-02:33:44 


1.1 (8.6) 


3.0 


0.27 (14.9) 


0.8 


0.051 (10.5) 


0.36 


T! 


SMM2 


18:13:37.6, 


-02:33:46 


0.8 (8.6) 


1.5 


0.14 (14.9) 


0.3 


0.033 (10.5) 


0.12 


V! 


CB 188 








8.4 




1.6 




0.30 


T5 


SMMl 


19:17:54.1, 


+ 11:30:02 


1.0 (9.3) 


2.0 


0.21 (15.8) 


0.3 


0.062 (10.5) 


0.10 


T5 


CB 199 


19:34:35.3, 


+07:27:24'' 


4.2 (8.6) 


15.1 


1.1 (14.9) 


4.0 








CB 205 








21.4 




3.5 




0.85 


T7 


SMMl 


19:43:22.1, 


+27:43:44 


2.4 (9.3) 




0.48 (15.8) 




0.078 (10.5) 




T7 


SMM2 


19:43:19.2, 


+27:43:19 


1.3 (9.3) 




0.30 (15.8) 




0.044(10.5) 




T7 


CB 224 








6.0 




0.8 




0.30 




SMMl 


20:35:30.6, 


+63:42:47 


1.2 (9.3) 


5.7 


0.27 (15.8) 


0.7 


0.067 (10.5) 


0.24 


■n? 


SMM2 


20:35:34.6, 


+63:43:10 


0.4 (9.3) 


0.4 


0.14 (15.8) 


0.1 


0.040(10.5) 


0.06 


TK 


CB 230 


21:16:53.7, 


+68:04:55 


3.9 (8.6) 


16.5 


0.92 (14.9) 


2.9 


0.23 (10.5) 


0.88 


TO 


CB 232 








27 




2.4 




0.45 


70 


SMM 


21:35:13.3, 


+43:07:13 


1.9 (12) 


3.5 


0.32 (14.9) 


0.9 


0.070(10.5) 


0.20 


■m 


IRSl 


21:35:14.5, 


+43:07:07 


1.6 (12) 


3.3 




0.3 




0.08 


70 


CB240 


22:31:45.9, 


+58:16:18 


<0.5' (8.6) 




0.14 (14.9) 


2.4 


0.028 (10.5) 


0.39 


TT 


CB243 












2.0 




0.68 


77 


SMMl 


23:22:52.1, 


+63:20:13 


1.6 (8.6) 


6.0 


0.37 (14.9) 


1.2 


0.083 (11) 


0.38 


72 


SMM2 


23:23:03.4, 


+63:20:20 


<0.5'^ (8.6) 




0.15 (14.9) 


0.5 


0.031 (11) 


0.19 


77 
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Table 5 — Continued 



Name 


Peak position" 


450;jm- 




850/jm 




1.3 mm 




Fig. 




(B1950) 






(fib) 


^loib 




^loib 








[ li : m : s, ° : ' : " ] 


[Jy/beam] 


[JyJ 


[Jy/beam] 


[Jy] 


[Jy/beam] 


[Jy] 






CB244 










5.1 




2.0 




23 




SMMl 


23:23:48.5, +74:01:08 


1.8 (8.6) 


13.0 


0.43 (14.9) 


1.6 


0.11 (10.5) 


0.8 




73 




SMM2 


23:23:29.5, +74:01:54 


<1.3' (8.6) 




0.32 (14.9) 


3.5 


0.06 (10.5) 


1.2 




73 




CB246 






30 




3.1 




2.56 




55 




SMMl 


23:54:11.0, +58:17:15 


1.3 (15) 


22 


0.12 (16) 


1.8 


0.039 (15) 


0.59 




25 




SMM2 


23:54:03.7, +58:18:36 


1.6 (15) 


8 


0.11 (16) 


1.3 


0.041 (15) 


0.58 































"Derived from gaussian fit to the component in the 1.3 mm map (if not stated otherwise). 

'^Flux density measured in a polygon enclosing the 2 cr contour. 

^3 cr detection limit; source not detected. 

*^ Source position derived from the 850/jm SCUBA map. 

'^The two sub-cores are not resolved in the 1.3 mm map. 



Table 6 

Globule morphology and multiplicity 



Name 


- Size" - 
['] [pc] 


Aspect 
Ratio 


A/ 
[mag] 


Globule 
Moiphology'^ 


Multiplicity of embedded sources 


Projected separations 
[AU] 


CB6 


5.5 


0.96 


2.6 


0.3 


1, T, N 


Unresolved 




CB 17 


2.2 


0.16 


1.0 


0.4 


I, T 


Two mm cores and IRS 


2,000 - 5,000 


CB26 


4.8 


0.20 


2.0 


1.5 


CO, C, BR 


Unresolved'^ 




CB34 


3.1 


1.35 


2.0 


1.0 


1, N 


Thi"ee mm cores and NIR star clusters 


6,000 - 80,000 


CB54 


4.4 


1.40 


1.6 


0.5 


1, N 


Unresolved mm core and NIR star cluster 


10,000 


CB58 


5.1 


2.22 


2.3 


1.0 


CO, BR 


Two mm cores and NIR star cluster 


7,500 - 50,000 


BHR12 


2.4 


0.28 


1.5 


1.2 


CO, BR, N 


Two mm cores 


8,000 


BHR36 


3.5 


0.41 


1.3 


1.5 


CO, BR, N 


Unresolved 




BHR55 


3.5 


0.30 


1.3 


1.2 


1, C 


Unresolved 




BHR71 


4.9 


0.25 


2.7 


3.4 


1 


Unresolved mm core, two IRS 


3,400 


BHR86 


7.7 


0.43 


1.7 


1.3 


CO, BR 


Unresolved 




CB 68 


3.9 


0.18 


1.3 


1.2 


1, T 


Unresolved 




BHR137 


3.9 


0.79 


1.7 


2.0 


I, N, T 


Unresolved mm core and IRS 


2,000 - 15,000 


CB 130 


3.5 


0.20 


2.5 


3.7 


I 


Two mm cores and IRS 


3,000 - 6,000 


CB 188 


1.6 


0.14 


2.0 


3.4 


I 


Unresolved 




CB 199 


5.0 


0.14 


1.2 


0.6 


I 


Unresolved 




CB205 


8.2 


4.76 


1.5 


1.5 


I, N 


Two mm cores and NIR star cluster 


4,000 - 60,000 


CB224 


3.4 


0.39 


1.0 


1.3 


I 


Two imn cores 


14,000 


CB230 


2.7 


0.31 


1.5 


1.5 


1, BR, N 


Unresolved mm core, two IRS 


4,000 


CB232 


3.5 


0.61 


2.5 


1.7 


CO 


MM core and IRS 


7,000 


CB 240 


2.7 


0.39 


1.5 


2.4 


I 


MM core and NIR star cluster 


1,200- 10,000 


CB243 


4.7 


0.96 


4.6 


2.0 


F, T 


Two mm cores and IRS 


7,000 - 50,000 


CB244 


7.9 


0.46 


2.0 


2.0 


I 


Two mm cores 


18,000 


CB246 


5.9 


0.24 


1.7 


2.7 


I, T 


Two mm cores 


13,000 



''Characteristic or mean angular size ( ^fcih), derived from the optical globule sizes listed by |Clemens & Barvainis||l988j for CB (noithern) 
globules and in Bourke et al. i 1995b i for BHR (southern) globules. Mean linear size derived by using the distances listed in Table|2] At the 
globule centres, the extinction is likely to be higher than this. 

'^Globule morphology, based on visual inspection of DSS2-red optical images: (1) isolated; (CG) cometary globule; (F) filamentary; (C) 
complex, multiple cores; (BR) bright rim; (N) nebula associated; (T) tail. 

^CB 26: contains possible a 20 AU binaiy with a circumbinary disk jUaunhardt et al.|2009|. 
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Fig. 1. — CB6: a) Optical image (DSS2 red) with contours of the 850 /jm dust continuum emission overlaid, b) NIR K-band image with 850 /im 
dust continuum contours. The IRAS PSC position is marked by a dashed ellipse, c) 450 jum dust continuum emission; grey scale and contours at 1.1 
and 1.5 Jy/beam. d) 850//m dust continuum emission; grey scale and contours at 35 to 175 by 35 mJy/beam. e) 1.3 mm dust continuum emission; 
grey scale and contours at 12 to 48 by 12mjy/beam. Beam sizes are indicated as grey ellipses. The cross marks the position of the 1.3 mm peak, 
f) Spectral energy distribution showing the IRAM 1.3 mm, SCUBA 850 and 450 jum, IRAS PSC, Spitzer MIPS and IRAC, and ground-based NIR 
data. 
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Fig. 2. — CB 17: a) Optical image (DSS2 red) with contours of the 1.3 mm dust continuum emission overlaid, b) NIR K-band image with 1.3 mm 
dust continuum contours. The IRAS PSC position is marked by a dashed ellipse, c) 450/jm dust continuum map; the source is not detected at rms 
400mjy/beam. d) 850/im dust continuum emission; grey scale and contours at 50 to 170 by 30mjy/beam. e) 1.3 mm dust continuum emission; 
grey scale and contours at 12 to 44 by 8mjy/beam. Beam sizes are indicated as grey ellipses. Crosses mark the positions of SMMl and SMM2 
in the 1.3 mm map. The asterisk marks the position of the Spitzer source IRS. f) Spectral energy distributions of SMM(l+2) and IRS, showing the 
IRAM 1.3 mm, SCUBA 850/<m, IRAS PSC, Spitzer MIPS and IRAC, and ground-based NIR data. Note that the 160 yum data points result only 
from a coarse flux splitting (85% IRS and 15% SMM) since the two sources are not resolved in the MIPS3 map. 
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Fig. 3. — CB26: a) Optical image (DSS2 red) with contours of the 850 //m dust continuum emission overlaid, b) NIR K-band image with 1.3 mm 
dust continuimi contours. The IRAS PSC position is marked by a dashed ellipse, c) 450 //m dust continuum emission; grey scale and contours at 
350, 700, 1050 to 3150 by 700mJy/beam. d) 850 //m dust continuum emission; grey scale and contours at 30, 60, 90, 150 to 400 by lOOmJy/beam. 
e) 1.3 mm dust continuum emission; grey scale and contours at 15, 30 to 150 by 30mjy/beam. Beam sizes are indicated as grey ellipses, f) Spectral 
energy distributions, showing the OVRO 3 mm, IRAM 1.3 mm, SMA 1.1 mm, SCUBA 850 and 450/<m, CSO 350 /im, IRAS PSC, Spitzer MIPS 
and IRAC, and ground-based NIR data. 
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Fig. 4. — CB 34: a) Optical image (DSS2 red) with contours of the 1.3 mm dust continuum emission overlaid, b) NIR K-band image with 1.3 mm 
dust continuum contours. The IRAS PSC position is marked by a dashed ellipse, c) 450 jmi dust continuum emission; grey scale and contovus at 500 
to 1700 by 300 mjy/beam. d) 850 fim dust continuum emission; grey scale and contours at 75 to 450 by 75 mjy/beam. e) 1 .3 mm dust continuum 

emission; grey scale and contours at 18 to 108 by 15 mJy/beam. Beam sizes are indicated as grey ellipses. Crosses mark the positions of the three 
1.3 ram peaks. No spectral energy distribution is shown since many different, partially unresolved sources contribute to the measured flux densities. 
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Fig. 5. — CB54: a) Optical image (DSS2 red) with contours of the 850//m dust continuum emission overlaid, b) NIR K-band image with 850/jm 
dust continuum contours. The IRAS PSC position is marked by a dashed ellipse, c) 450 jum dust continuum emission; grey scale and contours at 
0.4, 0.8, 1.2 to 9 0.8 Jy/beam. d) 850 /jm dust continuum emission; grey scale and contours at 50, 100, 200, 300, to 1500 by 200mjy/beam. e) 
1.3 mm dust continuum emission; grey scale and contours at 20, 40, 60, 90, 140, 190, 240, and 290mjy/beam. Beam sizes are indicated as grey 
ellipses. No spectral energy distribution is shown since many different, partially unresolved sources contribute to the measured flux densities. 
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Fig. 6. — CB58: a) Optical image (DSS2 red) witli contours of the 1.3 mm dust continuum emission overlaid, b) NIR K-band image with 
1.3 mm dust continuum contours. The IRAS PSC position is marked by a grey ellipse, d) 850^m dust continuum map; the source is not detected 
at rms 90mjy/beam. e) 1.3 mm dust continuum emission; grey scale and contours at 33 to 132 by 33 mjy/beam. Beam sizes are indicated as grey 
ellipses. No spectral energy distribution is shown since many different, partially unresolved sources contribute to the measured flux densities. 
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Fig. 7. — BHR 12: a) Optical image (DSS2 red) with contours of the 850/<m dust continuum emission overlaid, b) NIR K-band image with 
850 //m dust continuum contours. The IRAS PSC position is marked by a dashed ellipse, c) 450 jum dust continuum emission; grey scale and 
contours at 0.3, 0.6, 1 to 3.5 by 0.5 Jy/beam. d) 850 /jm dust continuum emission; grey scale and contours at 40, 100, 200 to 1200 by 150 mjy/beam. 
e) 1.3mm dust continuum emission; grey scale and contours at 25, 60, 100 to 600 by lOOmJy/beam. Beam sizes are indicated as grey ellipses. 
Crosses mark the positions of the two 850/jm peaks, t) Spectral energy distributions of SMMl (filled squares) and SMM2 (empty circles) showing 
the ATCA 3 mm, SEST 1 .3 mm, SCUBA 850 and 450 fim, IRAS PSC. Spitzer MIPS and IRAC, and ground-based NIR data. The 1.3 mm fluxes 
(sources not resolved in the SEST map) are derived by tilting two Gaussian sources to the map. 
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Fig. 8. — BHR 36: a) Optical image (DSS2 red) with contours of the 1.3 mmm dust continuum emission overlaid, b) NIR K-band image with 
1.3 mm dust continuum contours. The IRAS PSC position is marked by a dashed ellipse, d) 850jum dust continuum emission; grey scale and 
contours at 0.15, 0.3, 0.5, 0.75, and 1 Jy/beam. e) 1.3 mm dust continuum emission; grey scale and contours at 30 to 280 by 30mJy/beam. Beam 
sizes are indicated as grey ellipses, t) Spectral energy distribution showing the SEST 1.3 mm, SCUBA 850 /jm IRAS PSC, Spitzer MIPS and IRAC, 
and ground-based NIR data. 



36 




-12 



-14 






1 [ffln] 


1300 






1 ' 1 ' 1 ' 












.1 ■ 



IS 14 

la g f ^2 




30 -30 

AH, A. [oTEaBc] 




eo SO -SO 

Fig. 9. — BHR55: a) Optical image (DSS2 red) with contours of the 1.3 mmm dust continuum emission overlaid, b) NIR K-band image with 
1.3 mm dust continuum contours. The IRAS PSC position is marked by a dashed ellipse, e) 1.3 mm dust continuum emission; grey scale and 
contours at 55, 110, 165mJy/beam. The beam size is indicated as a grey ellipse, f) Spectral energy distribution showing the SEST 1.3 mm, IRAS 
PSC, Spitzer MIPS and IRAC, and ground-based NIR data. Note that the association of the IRAC fluxes with the FIR/mm source is somewhat 
uncertain. 
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Fig. 10. — BHR71: a) Optical image (DSS2 red) with contours of tlie 1.3 mmm dust continuum emission overlaid (see |Bourke et al. 11997} . 
b) NIR K-band image with 1.3 mm dust continuum contours. The IRAS PSC position is marked by a dashed ellipse, e) 1.3 mm dust continuum 
emission; grey scale and contours at 100, 250 to 1750 by 300mJy/beam. The beam size is indicated as a grey ellipse. Embedded infrared sources 
IRSl and IRS2 refer to Bourke (2001 1, t) Spectral energy distribution showing the ATCA 3mm, SEST 1.3mm, SMA 1.1 mm, SCUBA 850 and 
450/jm. CSO 350/im, IRAS 100 to 12^m, MIPS 70 and 24,um, IRAC 8 to 3.6/jm, as well as ISOCAM 15 and 7//m data. Note that this is a 
composite SED containing the sum of flux contributions from both IRSl and IRS2. 
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Fig. 1 1. — BHR86: a) Optical image (DSS2 red) witli contours of tlie 1.3 mmm dust continuum emission overlaid, b) NIR K-band image 
with 1.3 mm dust continuum contours. The IRAS PSC position is marked by a dashed ellipse, e) 1.3 mm dust continuum emission; grey scale and 
contours at 75 to 315 by 80mJy/beam. The beam size is indicated as a grey eUipse. f) Spectral energy distribution showing the SEST 1.3 mm, IRAS 
FSC, and Spitzer MIPS and IRAC data. 
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Fig. 12. — CB 68: a) Optical image (DSS2 red) with contours of the 850 //m dust continuum emission overlaid, b) NIR K-band image with 
850 /<m dust continuum contours. The IRAS PSC position is marked by a dashed ellipse, c) 450 /jm dust continuum emission; grey scale and 
contours at 1.2, 2.2, and 3.2Jy/heam. d) 850;um dust continuum emission; grey scale and contours at 90 to 540 by 90mJy/beam. e) 1.3mm dust 
continuum emission; grey scale and contours at 18 to 162 by 18 mjy/beam. Beam sizes are indicated as grey ellipses, f) Spectral energy distribution 
showing the OVRO 3 mm, IRAM 1.3 mm, SCUBA 850 and 450 jum, IRAS PSC and FSC, Spitzer MIPS and IRAC, and ground-based NIR data. 
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Fig. 13. — BHR 137: a) Optical image (DSS2 red) with contours of the 1.3 mmm dust continuum emission overlaid, b) NIR K-band image 
with 1 .3 mm dust continuum contours. The IRAS PSC position is marked by a dashed ellipse, e) 1 .3 mm dust continuum emission; grey scale and 
contours at 70 to 350 by 70 mjy/beam. The beam size is indicated as a grey ellipse. No spectral energy distribution is shown since it is not clear 
which of the sources emitting at shorter wavelengths, if any, is associated with the mm source. 
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Fig. 14. — CB 130: a) Optical image (DSS2 red) with contours of the 1.3 mm dust continuum emission overlaid, b) NIR K-band image with 
1.3 mm continuum contours, c) 450 jum dust continuum emission; grey scale and contours at 0.5 and I Jy/beam. d) 850 jum dust continuum emission; 
grey scale and contours at 45 to 270 by 45 mjy/beam. e) 1.3 mm dust continuum emission; grey scale and contours at 12 to 72 by 12 mjy/beam. 
Beam sizes are indicated as grey ellipses, f) Spectral energy distribution of SMMl (filled squares), SMM2 (submm only, empty squares), and IRS 
(NIR to MIR only, open circles), showing the IRAM 1 .3 mm, SCUBA 850 and 450 jxm, Spitzer MIPS and IRAC, and ground-based NIR data. 
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Fig. 15. — CB 188: a) Optical image (DSS2 red) with contours of the 850/jm dust continuum emission overlaid, b) NIR K-band image with 
850/<m dust continuum contours. The IRAS PSC position is marked by a dashed ellipse, c) 450/im dust continuum emission; grey scale and 
contours at 0.45 and 0.9Jy/beam. d) 850 jum dust continuum emission; grey scale and contours at 40 to 160 by 40mjy/beam. e) 1.3mm dust 
continuum emission; grey scale and contours at 15 to 60 by 15 mjy/beam. Beam sizes are indicated as grey ellipses, f) Spectral energy distribution 
showing the OVRO 3 mm, IRAM 1.3 mm, SCUBA 850 and 450 //m, IRAS PSC, Spitzer MIPS and IRAC, and ground-based NIR data. 
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Fig. 16. — CB 199 (B 335): a) Optical image (DSS2 red) with contours of the 850/im dust continuum emission overlaid, b) NIR K-band image 
with 850,um dust continuum contours. The IRAS PSC position is marked by a dashed ellipse, c) 450/im dust continuum emission; grey scale and 
contours at 0.4 to 4.4 by 0.5 Jy/beam. d) 850/im dust continuum emission; grey scale and contours at 60, 120, 200 to 1000 by 200mjy/beam. Beam 
sizes are indicated as grey ellipses, f) Spectral energy distribution showing the IRAM 1.3 mm jMotte & Andre|2001f , SCUBA 850 and 450/jm, 
IRAS PSC, and Spitzer MIPS and IRAC data. 
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Fig. 17. — CB 205: a) Optical image (DSS2 red) witli contours of tlie 850/im dust continuum emission overlaid, b) NIR K-band image witfi 
850 /<m dust continuum contours. The IRAS PSC position is marked by a dashed ellipse, c) 450/im dust continuum emission; grey scale and 
contours at 0.65 to 3.05 by 0.6 Jy /beam, d) 850/jm dust continuum emission; grey scale and contours at 120 to 470 by 70mJy/beam. e) 1.3 mm dust 
continuum emission; grey scale and contours at 16 to 64 by 16mjy/beam. Beam sizes are indicated as grey ellipses. Crosses mark the positions of 
the two main 1.3 mm peaks. No spectral energy distribution is shown since many different, partially unresolved sources contribute to the measured 
flux densities. 
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Fig. 18. — CB 224: a) Optical image (DSS2 red) with contours of the 1.3 mm dust continuum emission overlaid, b) NIR K-band image with 
1.3 mm dust continuimi contours. The IRAS PSC position is marked by a dashed ellipse, c) 450 /<m dust continuum emission; grey scale and 
contours at 0.24 to 1 .44 by 0.24 Jy/beam. d) 850 yum dust continuum emission; grey scale and contours at 60 to 240 by 60 mjy/beam. e) 1 .3 mm dust 
continuum emission; grey scale and contours at 12 to 60 by 12 mJy/beam. Beam sizes are indicated as grey ellipses. Crosses mark the positions of 
the two 1 .3 mm peaks, f) Spectral energy distribution of SMMl (filled squares) and SMM2 (open circles) showing the OVRO 3 mm, IRAM 1.3 mm, 
SCUBA 850 and 450 ^/m, IRAS PSC and FSC, and ground-based NIR data. 
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Fig. 19. — CB 230: a) Optical image (DSS2 red) with contours of the 1.3 mm dust continuum emission overlaid, b) NIR K-band image with 
1.3 mm dust continuum contours. The IRAS PSC position is marked by a dashed ellipse, c) 450 ;<m dust continuum emission; grey scale and 
contours at 350, 700, 1050 to 3850 by 700mJy/beam. d) 850/im dust continuum emission; grey scale and contours at 30, 60, 90, 150 to 900 by 
lOOmJy/beam. e) 1.3 mm dust continuum emission; grey scale and contours at 15, 30 to 210 by 30mjy/beam. Beam sizes are indicated as grey 
ellipses, f) Spectral energy distributions of IRSl (filled squares) and IRS2 (open circles, NIR to MIR only), showing the OVRO 3 mm, IRAM 
1.3 mm, SCUBA 850 and 450yum, UKT14 350 fim, IRAS PSC and FSC, Spitzer MIPS and IRAC 8 to 3.6//m, ISOCAM 6.7 /jm, and ground-based 
NIR data. Note that contributions from both IRSl and IRS2 are included in the submm SED of SMMl because the source is not resolved at 
wavelengths longward of 10/zm. However, IRSl dominates at all wavelengths. 
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Fig. 20. — CB 232: a) Optical image (DSS2 red) with contours of the 850/im dust continuum emission overlaid, b) NIR K-band image with 
850 fim dust continuum contours. The IRAS PSC position is marked by a dashed ellipse, c) 450 pm dust continuum emission; grey scale and 
contour at 1.4mJy/beam. d) 850/im dust continuum emission; grey scale and contours at 70 to 320 by 50mJy/beam. e) 1.3 mm dust continuum 
emission; grey scale and contours at 12 to 60 by 12mjy/beam. Beam sizes are indicated as grey ellipses. The cross marks the position of SMMl 
in the 1.3 mm map. Asterisks mark the position of the two nearby Spitzer IRS. f) Spectral energy distribution of CB 232 IRSl (filled squares) and 
SMM (open circles, submm only), showing the IRAM 1.3 mm, SCUBA 850 and 450 jum, IRAS PSC, Spitzer MIPS and IRAC, and ground-based 
NIR data. 
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Fig. 21. — CB 240: a) Optical image (DSS2 red) with contours of the 1.3 mm dust continuum emission overlaid, b) NIR K-band image with 
1.3 mm dust continuum contours. The IRAS PSC position is marked by a dashed ellipse, c) 450 /im dust continuum map; the source is not detected 
at rms 170mjy/beam. d) 850 /im dust continuum emission; grey scale and contour at lOOmJy/beam. e) 1.3 mm dust continuum emission; grey 
scale and contours at 16 and 32 mJy^eam. Beam sizes are indicated as grey ellipses. No spectral energy distribution is shown since many different, 
partially unresolved sources contribute to the measured flux densities. 
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Fig. 22. — CB 243: a) Optical image (DSS2 red) with contours of the 1.3 mm dust continuum emission overlaid, b) NIR K-band image with 
1.3 mm dust continuum contours. The IRAS PSC position is marked by a dashed ellipse, c) 450 /jm dust continuum emission; grey scale and 
contours at 0.5 to 2 by 0.5 Jy/beam. d) 850 ;um dust continuum emission; grey scale and contours at 60 to 360 by 60mJy/beam. e) 1.3 mm dust 
continuum emission; grey scale and contours at 15 to 75 by 15mjy/beam. Beam sizes are indicated as grey ellipses. Crosses mark the positions 
of of SMMl and SMM2. The asterisk marks the position of the reddened star and IRAS point source, f) Combined spectral energy distribution of 
SMMI and IRS, showmg the IRAM 1.3 mm, SCUBA 850 and 450jt(m, IRAS PSC, Spitzer MIPS, MSX, and ground-based NIR data. Fluxes are 
dominated by SMMl at /I > 450 /<m and by IRS at /I < 100 jum, but the two SEDs could not be separated. 
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Fig. 23. — CB 244: a) Optical image (DSS2 red) with contours of the 1.3 mm dust continuum emission overlaid, b) NIR K-band image with 
1.3 mm dust continuimi contours. The IRAS PSC position is marked by a dashed ellipse, c) 450 /<m dust continuum emission; grey scale and 
contours at 0.6 to 2.2 by 0.5 Jy/beam. d) 850 /<m dust continuum emission; grey scale and contours at 100 to 500 by 80mjy/beam. e) 1.3 mm dust 
continuum emission; grey scale and contours at 15 to 105 by 15 mjy/beam. Beam sizes are indicated as grey ellipses. Crosses mark the positions 
of SMMl and SMM2. f) Spectral energy distribution of CB 244 MM 1 showing the OVRO 3 mm, IRAM 1.3 mm, SMA 1.1 mm, SCUBA 850 and 
450 //m, UKT14 350 ^^m, IRAS PSC and FSC, Spitzer MIPS and IRAC, and ground-based NIR data. 
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Fig. 24. — CB246: a) Optical image (DSS2 red) witii contours of the 1.3 mm dust continuum emission overlaid, b) NIR K-band image with 
1.3 mm dust continuum contours, c) 450;um dust continuum emission; grey scale and contour at 1.0 Jy/beam. d) 850,i;m dust continuum emission; 
grey scale and contours at 60 and lOOmJy/beam. e) 1.3 mm dust continuum emission; grey scale and contours at 15, 25, and 35 mjy/beam. Beam 
sizes are indicated as grey ellipses. Crosses mark the approximate central positions of the two main submm cores. No spectral energy distribution 
is shown because no fluxes are detected at wavelengths shortward of 450 //m. 
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Fig. 25. — Compilation of SEDs (same as in Figs .[l]tlirough [23) , ordered by Tboi and evolutionary stage (of. Table^. Note that the classification 
of the last two sources (CB 17-IRS and CB 130-IRS) is uncertain because of missing long-wavelengths constraints to the SED. 
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Fig. 26. — Ratio of submm to bolometric luminosity vs. bolometric temperature for the sources compiled in Table|7] The vertical dashed line 
marks the Tbol = 70 K boundary between Class and Class I sources proposed by [Myers et al.| Jl998j. Sources with NIR/MIR nebulosity only 
(i.e. no star-like point source) are marked as filled squares. Sources with NIR point source (star) are marked as open asterisks. See Sect. |4. l| for a 
discussion of evolutionary tracers. 
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Fig. 27. — Mean I^WHM core size and source-averaged density in globule cores as function of evolutionary stage (PS = prestellar core. Class 
protostars. Class 1 YSO cores and extended envelopes). Note that most Class I soitfces exhibit unresolved (sub)mm sources with presumably 
significant flux contributions from circumstellar disks. 
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Fig. 28. — Four examples of globule cores in which MIR and submm sources are not identical and which obviously contain two or more 
sources of different evolutionary stage with a few thousand AU: Top left: CB 17, IRACfiiim image (greyscale) and 1.3 mm dust continuum contours 
(see Fig.|2). Top right: CB 130, IRAC 8 fim image (greyscale) and 1.3 mm dust continuum contours (see Fig. |14) . Bottom left: CB232, MIPS 
24 /im image (greyscale) and 850/im dust continuum contours (see Fig. |20( . Bottom right: CB243, MIPS 24 /im image (greyscale) and 850//m 
dust continuum contours (see Fig. |22|. I RAS PSC poitions are marked by dashed error ellipses, where applicable. Linear scale bars refere to the 
distances listed in Table|2] See Sect. l4.2l for discussion. 
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